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Introduction 



The success of Feynman's parton model [1] has pushed the search of a field the- 
ory formulation of the parton hypothesis. The discovery of asymptotic freedom 
of unbroken non-abelian gauge theories, which implies that the corresponding 
particle interaction weakens at short distances, has clearly indicated the held 
content of the Feynman model. Partons are fermions, quarks, and vector bosons 
interacting through exchanges of vector bosons, gluons. 

Asymptotic freedom [2] makes more precise Feynman's idea that, whenever 
one has high energy and momentum transfer, partons behave as free particles. 
This is indeed what happens in Deep Inelastic Scattering (DIS) and in quark 
anti-quark annihilations into heavy bosons. 

However, partons, therefore quarks and gluons, should be bound in a 
hadronic initial state. Therefore the partonic process should be characterized 
by two scales, the high energy momentum transfer scale and the binding energy 
scale, which should be less than few hundred McV. 

Therefore if one tries to exploit perturbative field theory in order to compute 
corrections to the naive free parton model, one finds typically logarithmic factors 
in the ratio of the two scales, which is about one thousand. These logarithmic 
factors appear as multipliers of the strong structure constant a s . Therefore the 
validity of perturbation theory becomes less obvious in much the same way as 
the idea of parton independence, that is, factorization. Feynman amplitudes 
in these calculations must be computed in the neighbourhood of singularities 
which are typically of the collinear type. The only simplification induced by the 
presence of the two above mentioned scales is that at least for the quarks up and 
down, whose mass is few MeV, one can consider quarks as massless as gluons are. 
Therefore the field model should be considered in first (free) approximation scale 
invariant and the radiative corrections should violate Bjorken scale invariance 
due to renormalization. This is indeed what happens and the corrections to 
Bjorken scaling [3] give a strong evidence in favour of the non-abelian QCD 
model. 

However the calculations of radiative corrections together with the presence 
of the above mentioned scales, one of which is associated with the process, the 
other appearing in the hadronic wave function, have posed non-trivial technical 
difficulties which have been overcome thanks to the discovery of dimensional 
rcgularization [4, 5] and of the corresponding minimal renormalization schemes. 

The simplicity of computation is only one of the advantages which have 
strongly favoured the use of the dimensional method. A second very important 
advantage is that it doesn't break gauge invariance and applies without problems 
to the quark and gluon mass-shell amplitudes. We shall see that this is very 
important in view of parton factorization hypothesis. The method is consistent 
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with the description of the hadron as a gas of free quarks and gluons whose 
distribution function has to be computed on the basis of consistency conditions 
of different scale choices. 

Even if the free gas scheme is not physically consistent, otherwise the 
hadronic density of state should be exceedingly large, the presence of a consis- 
tent computational method for high energy processes has given solid technical 
foundations to the parton model. 

In this framework a huge amount of work has been done aiming at the 
analysis of the singularities of Feynman diagrams as functions of the kinematic 
invariants, that is, of the invariant squared partial sum of the momenta carried 
by the external vertices. In particular, in the study of DIS, the behaviour of a 
mass-shell QCD amplitude when partons are considered massless is related to 
collinear singularities of the amplitude. The study of these singularities gives 
information through the Altarelli-Parisi approach [6] on the parton distribution 
function of e.g. a target baryon. 

Collinear singularities appear in dimensional regularization of a mass-shell 
massless parton amplitude when the dimensional regularization parameter e = 
goes to zero; that is collinear singularities appear as poles in the origin 
of the e complex plane in much the same way as UV singularities do. This 
implies that one has to avoid IR-UV singularity mixings by first subtracting 
UV divergences. 

An important result concerning the Feynman graph singularities is that in 
the Schwinger parametric representation [7] the divergent parts of a diagram 
are confined in particular sectors of the parametric space. 

In fact in a mass-shell amplitude, using dimensional regularization, one finds 
that e — >• singularities associated with collinear divergences appear only in 
the restriction of the parametric representation to particular sectors which can 
be identified using the results published by Speer in 1975 [8]. This induces a 
remarkable simplification in the singularity analysis but until recently it has not 
been exploited in calculations. (See however [9] and references therein). 

On account of these results and of more recent analogous formulations [10], 
it is reasonable to verify the simplifying power of Speer's results in the calcu- 
lation of collinear singularities of parton processes. However, beyond purely 
technical aspects, an easy tool for the computation of the infrared singularities 
of the Feynman amplitudes allows a better physical insight into Parton Physics. 
From the point of view of the physical interpretation, one should take into ac- 
count that partons are identified with bound, and hence off-shell quarks and 
gluons. However single parton amplitudes with off-shell initial particle states 
present a further technical difficulty. Indeed it is fairly well known since the 
seminal Yennie-Frautschi-Suura work [11] that, already in QED, off-shell elec- 
tron amplitudes are not "gauge invariant" and that this has consequences on 
the compensation of infrared (soft) divergences of the single graphs. Here gauge 
invariance should be given two different meanings, the most common one being 
the vanishing of space-time divergence of any conserved current vertex; there is 
however a second meaning concerning the independence of the gauge fixing of 
quantization choice. As a matter of fact strictly speaking the single parton con- 
tributions are not physically meaningful and hence parton factorization, which 
is the basis of the parton model, should be taken " cum grano salis" . 

As a matter of fact hadrons are colour neutral and hence the cross section of 
a fully inclusive deep inelastic process is free of infrared singularities. Indeed this 
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is a consequence of the famous Kinoshita-Lee-Nauenberg theorem [12], which 
can be easily understood noticing that the cross section is the absorptive part of 
a virtual forward Compton scattering whose amplitude should not be affected 
by infrared divergences. 

There are however two obstructions in pushing forward this point of view. 
The first one is that this would bring the analysis far away from the independent 
parton model, the second one is the need of knowledge of the hadron-parton 
vertex which cannot be deduced from perturbative theory. 

One could try to overcome these difficulties assuming that the partons be 
" almost" on-shell and the collinear singularities be induced by the " almost" on- 
shellness. However the analysis should also consider the singularities due to the 
exact on-shell nature of the final state partons. 

In QED, when the initial state is neutral, final state singularities come from 
the coupling of real and virtual soft photons to the final charged particles and 
cancel in the inclusive cross sections. This happens also in first order corrections 
to QCD amplitudes and suggests that the final parton on-shell singularities 
should be forgotten keeping only the contributions which are singular when the 
initial parton goes on-shell. 

Notice that considering off-shell amplitudes is a kind of regularization since 
the infrared singularities are regularized by a suitable off-shell (Euclidean) choice 
of external momenta and are therefore well separated from the (dimensionally 
regularized) UV singularities. Indeed studying a massless theory for space-like 
external momenta [8, 13], one shows that if the external momenta are non- 
exceptional, i.e. no partial sum vanishes, the diagram has no infrared diver- 
gences if the field have positive mass dimension and all the internal vertices 
have dimension four, as in QCD. 

Hence the idea is to regularize the parton amplitudes choosing off-shell initial 
momenta [14] and to select the contributions which are singular when the parton 
goes on-shell. Notice that contrary to the off-shell amplitude, its mass-shell 
singular part should be gauge invariant and independent. 

We hope that applying the present idea to a physical situation, such as Deep 
Inelastic Scattering, will help us to make this point clear. 

In this work we consider the difficult off-shell analysis in the light of few, 
more or less, recent progresses in Feynman graph computations. These are 
essentially based on the extended use of Mellin-Barnes transform [15], aiming 
at the singularity analysis, and the sector decomposition of parametric Feynman 
integrals, mentioned above. 

It is clear that a generic not renormalized amplitude could present UV sin- 
gularities which could combine with mass singularities. The presence of UV and 
mass singular contributions could appear in the case of multiloop diagrams with 
UV divergent subdiagrams. However, once these ultraviolet divergences are sub- 
tracted using the Breitenlohner-Maison scheme [16], they do not interfere with 
the mass singularities [17]. 

The best way of describing a computational algorithm consists in applying 
it to a simple, however not trivial case. What we are going to show in this thesis 
is how the combined sectorialization-Mcllin Barnes transform technique works 
in the study of collinear singularities in the first order in the strong constant a s 
structure functions of DIS. 
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Chapter 1 

Massless QCD 



Currently Particle Physics is successfully described by the Standard Model [18, 
19], i.e. the theory of strong and electroweak interactions, based on gauge group 
SU(3) ® SU(2) ® U(1)y- The corresponding gauge fields are gluons, which are 
the mediators of the strong interaction, and the electroweak bosons, the photon, 
W ± and Z . 

The Standard Model contains Quantum Chromodynamics or QCD, the non- 
abelian gauge theory which describes strong interactions based on the colour 
group SU(3). 

SU(2)®U(1)y is the Glashow-Weinberg-Salam electroweak symmetry group. 
The U(l)y factor is the abelian phase group generated by the weak hyperchage 
operator Y , which is related to the electric charge Q and to the weak isospin 
third component T 3; i.e. one of the SU{2) generators, by means of Q = T 3 + Y. 

W^ 1 and Z bosons, differently from the photon and the gluon, are massive; 
Lagrangian mass terms are induced by the partial spontaneous breakdown of the 
symmetry group SU{2) x U{1)y according to the Higgs mechanism. 

The unbroken gauge symmetry group is SU{2>) ® U(1)q. It corresponds to 
the exact symmetry of Nature; U(1)q is the phase group generated by the elec- 
tric charge. Its gauge boson is the photon. SU(3) gauge bosons are gluons: 
the eight gluons mediate the strong interactions between quarks, antiquarks and 
gluons, commonly called "the partons" . As a matter of fact quarks are massive, 
however, in particular up and down quarks have masses of few Me V. If one 
renormalizes the theory at the light quark mass scale pL in the physical ampli- 



tudes the coupling ct s appears multiplied by log ( — ) , where E corresponds to 



the energy scale of the process. If e.g. E is about 100 GeV, log(^) « 12 and 
hence the perturbative expansion parameter is too large. If, on the contrary, one 
chooses renormalization scale corresponding to E, the universally subtracted off- 
shell amplitudes become quark mass independent and hence, in general, they can 
be identified with the massless quark amplitudes 

In presence of massless particles, Feynman field theory, i.e. diagrammatic 
perturbative theory, produces infrared divergences, which necessarily have a phys- 
ical nature: they must be carefully studied in the light of the physical process. 

In conclusion: fundamental physics must deal with infrared divergence prob- 
lem because of the existence of massless particles. Infrared divergences arise 
particularly in QCD. In this chapter we give a brief summary of perturbative 
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QCD, focusing on applications in which infrared divergences appear; in this way 
we want to give the physical background of our work. 

1.1 Introduction to QCD 

QCD or Quantum Chromodynamics [20, 21, 22] is a gauge theory associated 
with the colour group SU(3), i.e. the group of 3 x 3 complex unitary matrices 
for which the determinant equals unity; they satisfy 

U+U = l, det(U) = l. (1.1) 

QCD is the mathematical model which describes strong interactions, that are 
generated by exchanges of gluons and quark-antiquark pairs in the hadrons. 
The SU (3) group acts on the fields of the theory, quarks, which correspond to 
colour triplets and gluons, corresponding to octets. 
The corresponding Lagrangian density is: 

C = -\i F aK + f^iiiPa - ™/)V>/, (1.2) 

A=l f=l 

where: 

• ipj is the spinor field associated with the j colour component of a quark 
with flavour / and mass m/. 

• p = -0^7^, where 7^ are Dirac matrices and D p is the covariant derivative: 

Dfx^j — d^Sij ie s A^ij. (1-3) 

In this last expression the gauge constant e s has been introduced. In 
general one refers to the strong coupling constant a s = the QCD 
constant in analogy with the fine structure constant in QED. 

One has introduced also A Mji j = ^2 A tAij • A A where A = [1,8]: A^ are 
the eight gluonic massless fields and t A are the eight SU(3) generators of 
the quark tridimensional representation; therefore they are a basis of the 
linear space of the 3x3 traccless matrices acting on quark ipj and they 
satisfy Tr (t^s) = Sab/2. Structure constants /abc are defined by the 
commutation relations [f^ts] = ifABctc- 

• F A V is an antisymmetric tensor that transforms under the gauge adjoint 
representation 

PV = V 1 with F„„ = (1.4) 

A 

and in components: 

F% - dpA? d v A$ e s f ABC A*AC. (1.5) 

It is believed that the very simple Lagrangian (1.2) describes the whole strong 
interaction physics: from the hadronic spectrum to confinement properties, to 
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high energy processes. In ordinary situations quarks can't be free, because of 
confinement, i.e. the existence of a quark-antiquark potential that rises lin- 
early with distance; free particles correspond only to states that are invariant 
by colour transformations. On the contrary high energy phenomena are char- 
acterized by asymptotic freedom: at short distances (and at high energies) the 
renormalized coupling constant becomes small, in contrast to what one observes 
in QED: therefore quarks inside hadrons appear nearly free. 

Asymptotic freedom is a property of utmost importance because it allows 
perturbative calculation in QCD processes characterized by short distances or by 
large energy or momentum transfer, such as "hard" processes. The perturbative 
method is very predictive as far these processes are concerned, particularly in 
electron-positron high energy annihilation into quark-antiquark pairs and in 
deep inelastic scattering. In these processes the comparison between theory and 
experiments provides a remarkable confirmation of QCD predictions. 

There is another method which allows the study of QCD in a non- 
perturbative way: this is lattice QCD, i.e. QCD on a finite and discrete space- 
time. Lattice calculations allow to compute physical observables by means of 
numerical simulations and to study QCD also at large coupling constant. In 
our work we consider only the perturbative approach, because it is the most 
used analytical method in experimental data analysis in the physical conditions 
when Hi-singularities become important. 

Starting from Lagrangian density given in equation (1.2), which describes 
classical QCD, one builds up Feynman rules corresponding to diagram vertices 
and propagators. QCD Feynman diagrams have different kinds of vertices since 
QCD is a non-abelian theory: in addition to the gluon-quark-antiquark vertex 
(similar to the QED photon-electron-positron vertex) one has also three and four 
gluons vertices and ghost field vertices; the ghost fields, which are unphysical 
scalar fields, are introduced into the theory in the Faddeev-Popov quantization 
procedure. 

In figure (1.1) QCD Feynman rules are shown. 

Using Feynman rules, one builds up the perturbative expansion of transi- 
tion amplitudes. One knows that Feynman diagrams with closed meshes, loops, 
correspond to integrals on the momentum space; such integrals could produce 
ultraviolet divergences in the high momentum region. In presence of massless 
particles and in the small momentum region, also infrared divergences appear. 
A regularization process is introduced to remove ultraviolet divergences pre- 
serving gauge symmetry and Lorentz invariance; this process yields UV-finite 
amplitudes. 

In general one computes cross sections which are related to the massless 
square of amplitudes which however are in general linear combinations of diffi- 
cult terms. More frequently one is interested in inclusive cross sections which 
can be obtained from the imaginary (absorptive) part of forward amplitudes. 
These are typically, when the final state is a many particle state, loop ampli- 
tudes, that is, amplitudes corresponding to diagrams with loops. 

Inclusive cross sections are in general less singular in their dependence of 
kinematic parameters than exclusive ones, those associated with a restricted 
selection of possible final states, since many potentially singular contributions 
compensate each other. 

There are many ways of proceeding. The main one is the dimensional regu- 
larization [4, 5]. 
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Figure 1.1: QCD Feynman rules. From the top to the bottom one has: gluon propaga- 
tor (expressed in Feynman gauge), ghost propagator and quark propagator. Then one 
has three gluon vertex, four gluon vertex, non-physical gluon-ghost-antighost vertex 
and gluon-quark-antiquark vertex. 
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Integrals are made finite varying space-time dimension from d = 4 to d = 
4 — 2e, where e is infinitesimal. Feynman amplitudes become meromorphic 
functions of the variable d and therefore of e. So ultraviolet divergences appear 
in the amplitude as poles in e = 0. Finite parts are obtained by subtracting 
Maclaurin expansion terms in e. 

Formally this subtraction is obtained by means of field and parameter re- 
definition (for example of e s ) which depend on e and have poles that cancel 
diagram divergences. This process is called renormalization. 

The quantum theory is represented by a Lagrangian density, called renor- 
malized, which has the same structure as the classical one. The redefinition of 
the starting Lagrangian parameters, called bare, in terms of physical fields and 
parameters and of e, allows to obtain amplitudes at every order of perturbation 
theory without ultraviolet divergences and to express amplitudes by means of 
physical parameters measured in experiments. 

1.2 Massless QCD and infrared divergences 

Equation (1.2), as we said, is the QCD classical Lagrangian density. Quark 
masses are the only parameters with physical dimensions, in particular with 
energy dimension using the natural system K = c = 1 . 

Therefore one expects that massless QCD is scale invariant. Scale invariance 
is the property for which observables are not dependent on the energy scale, but 
only depend on ratios of dimensional quantities. 

In reality QCD is scale invariant only classically and not after quantization 
because scale symmetry is broken by regularization and renormalization pro- 
cesses. Going from 4 dimensions to 4 — 2e dimensions, it is necessary to intro- 
duce a new dimensional parameter in the quantum version of the theory which 
introduces a new mass scale into the theory. In other words scale invariance is 
broken by regularization. Anyway deviations from scaling are asymptotically 
small, logarithmic and computable. 

In massive QCD there are other scaling corrections due to the presence of 
quark masses. However mass corrections decrease as powers of the ratios be- 
tween the masses and the energy scale in high energy processes and in particular 
in hard processes. 

Hard processes are characterized by energy and momentum exchanges much 
larger than Aqcdj the characteristic scale introduced by the breakdown of the 
scale invariance: 



As we said, hard process amplitudes, when masses are null or negligible, are 
affected by infrared divergences (IR) . These must be distinguished between soft 
infrared divergences and collinear ones, also called mass singularities. 

The first ones arise when the four-momentum of a real or virtual particle 
become soft, i.e. its components tend to zero. Finally one has mass or collinear 
singularities when one has transitions from a state of a massless particle to a 
state of two massless particles with four-momentum p^ and p ,f * moving parallel 
to each other. The invariant mass of the system becomes zero, although none 
of the momenta arc soft: 



Ei ~ E, E » Kqcd » mi . 



(1.6) 



k 2 = (p + p'f = 2EE'{\ - cos 6)->0 



(1.7) 
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k 



p+k p 

Figure 1.2: Diagram containing both a soft divergence and a mass singularity. 

In literature one defines collinear the divergences due to a term ( 1 _ c 1 os6) ) 
inside phase space integral (in cross section calculation); if instead the divergence 
appears inside a loop integral, it is called mass singularity. In this work we use 
two terms refer for the same thing. 

For example we consider a virtual quark (figure 1.2) that transforms into 
a real quark emitting a gluon at an angle 6 with respect to its momentum. 
The virtual quark propagator, which carries four-momentum p + k and which 
corresponds to the amplitude of the diagram in figure, is: 

1>+$ i (lg) 



(p + k) 2 -m 2 2(p-k) 2E k E p l-/3 p -cos0' 

where E k corresponds to the photon energy, E p and m are respectively quark 

energy and mass and one has j3 p = ^Jl — Since the gluon is massless, E k 

could vanish giving a soft infrared divergence. But when quark mass is zero, 
one has a second type of divergence, the collinear singularity. Indeed for m — > 
one has j3 p —¥ 1, therefore (1 — j3 p ■ cos 9) is null in the forward direction. Due to 
the presence of an helicity constraint which forbids the transition in the forward 
directions since angular momentum conservation requires a fermion helicity flip 
which is forbidden by the collinear helicity conservation, the transition proba- 
bility turns out to be proportional to 1 _l os0 and hence its phase space integral 
of the cross section, diverges as 



/ 



dcosd HQ) 
1-COS6T 1 ' 



This is a collinear divergence. 

A renormalization theory, which is able to remove infrared divergences, as in 
the case of ultraviolet singularities, doesn't exist. For example if one assigned 
a mass to the gluon, such that the amplitude would converge, one could not 
remove the mass dependence from the amplitude redefining the Lagrangian 
parameter. 

Nevertheless there are two theorems [12, 23] which in certain conditions 
guarantee the cancellation of the infrared divergences in physical processes. 

This suggests the idea that infrared divergences are due to an improper 
definition of initial and final states of the process. 

The first theorem is Bloch-Nordsieck's theorem, which is valid in QED and 
asserts that summing up over the final states in a little energy interval, soft in- 
frared divergences disappear because there is a cancellation between the virtual 
diagram divergences and the real ones. One understands this deletion, observ- 
ing that a particular scattering process, which consists in a transition from an 
initial state \a > to a final state \b >, is experimentally indistinguishable from 
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Figure 1.3: Electron-positron annihilation with hadrons production by means of an 
exchanged virtual gauge boson. 



the one with an arbitrary number of photons in the final state, provided that 
the photon energy is small enough. 

In QCD, this theorem is violated, but this happens because of subdominant 
terms, as it has been shown in [24] and in [25]. 

The second theorem is the Kino shita- Lee- Nauenberg (KLN^s one, which is 
valid in perturbation theory: collinear divergences, originating from massless 
particles in initial and/or final states, cancel up if one sums the transition prob- 
abilities between all of the states that are approximately mass degenerate with 
respect to the considered ones. 

Whereas the sum on degenerate final states is motivated by the finite res- 
olution of measure instruments and is an example of the already mentioned 
cancellation of singularities in inclusive cross sections, it is more difficult to ex- 
plain the sum on initial states: indeed one must first define the distribution of 
degenerate states upon which the sum must be carried. In QCD the condition 
which allows the cancellation of mass singularities coming from initial states 
generates Altarelli-Parisi equations. Mass singularities due to initial states are 
absorbed into the initial partonic densities, that correspond to the probability 
densities of finding the given parton into the initial hadron. Partons are the 
fundamental constituents, i.e. quarks and gluons. 

Considering hard processes in a more detailed way, one describes as an ex- 
ample the process of hadrons production coming out from electron-position 
annihilation. The reference is figure (1.3). 



cr(e + e — > hadrons) 
o- p oint(e + e- fi+IJ,-)' 



Historically the measure of this observable has confirmed the number of QCD 
colours, which is equal to 3, as we know. In (1.10) the point cross section is 

given by cr point (e + e~ — > = where s = Q 2 = AE 2 is the square of 

the center of mass (the invariant square of the total 4-momentum). 



13 



1.3 An example of hard process with infrared 
singularities: the Deep Inelastic Scattering 
(DIS) 



Figure 1.4: Deep Inelastic Scattering Diagram. 

Another important class of hard processes is the one known as Deep Inelastic 
Scattering (DIS) [21, 26], which groups processes of the type: 

l + N^-l' + X, l = e ± , /x ± , v, v. (1.11) 

DIS processes have a great importance in physics history because they clearly 
confirmed QCD predictions on the structure of nucleons, which are considered 
composite states of elementary particles: quarks and gluons. 1 

In figure (1.4) the DIS Feynman diagram is shown. 

Introducing the four-momentum transferred by the electron 

q = k — k', 

one defines kinematic invariants v = ^2, and Q 2 = —q 2 . Then one defines the 
Bjorken variable: 

O 2 

x b = f-, (1.12) 

which is adimensional and whose values belong to the interval between and 
1. In the case of elastic scattering one has Xb = 1; that is, Xb = 1 corresponds 
to a single parton final state. In the laboratory system, where the nucleon with 
mass m is at rest, lepton and nucleon four-momenta are: 

fc p = (E, k) ^ = (E', k>) p M = (m, 0, 0, 0), 

therefore the energy transferred from the electron to the nucleon is 

v = E-E'. (1.13) 

In covariant form the Bjorken variable is: 



1 The parton model treats hadrons as particles made of elementary constituents, said par- 
tons. Together with the study of the hadronic spectrum and of QCD itself, one assumes that 
partons are just quarks and gluons. 
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In case of deep inelastic scattering, the process duration, as it is seen by partons, 
is of order (\jv) and it is much less than the interaction time of partons, whose 
order of magnitude is the inverse nucleon mass: 

m « v. (1-15) 

In such conditions one can neglect quark interactions. 

In DIS experiments physicists measure gluon densities and quark densities 
as function of xb and of ^3-, where /i is the reference energy, they measure also 

a s (^j-), i.e. how the coupling constant depends on the momentum transfer. 

The cross section of the reaction is given by the product of a leptonic tensor, 
easy to evaluate because the lepton is point-like, and a hadronic tensor: 

a = L ati W afi . 

The hadronic tensor W a p is given by the Fourier transform of the invariant 
expectation value of the product of the two electroweak currents in the suitably 
normalized nucleon state: 



i(p,q)~ J dx cxp(iqx) < p | J+(x)J u (0)\p > . (1-16) 



This expectation value is directly related to the total (inclusive) cross section 
of the interaction of a nucleon with a virtual photon and hence, due to the optical 
theorem, to the imaginary part of the forward virtual photon- nucleon Compton 
scattering. 

In the case of an unpolarized nucleon current conservation and CP invariance 
imply that this tensor has the form: 

W aP (p,q) = -(g a( , + ^) Fl (x B ,Q 2 ) + ( Pa + Js-)(p p + JL. ) f »(xb,Q?) _ (117) 
Q z 2x B 2x B p ■ q 

Here the coefficient functions Fi(xB,Q 2 ) are called electromagnetic structure 
functions of the nucleon and are related to the hadronic tensor in the following 
way: 

2x B x B 2x B 

p a p* 3 W a , l3 (p,q) = ^(F2(x B ,Q 2 )-2x B F 1 (x B ,Q 2 )) = ^F L (x B ,Q 2 ),(l.l8) 

where Fl(xb,Q 2 ) is known in literature as longitudinal structure function. 

In the limit of very high energies Q 2 >> m 2 , i.e. in the limit in which the 
hadronic masses are negligible, the structure functions approximatively obey 
Bjorken scaling: 

F^xbM 2 ) -tFiixg) (1.19) 
F 2 (x B ,Q 2 ) -^F 2 (x B ). 

In other words in the above limit, i.e. for large values of Q 2 , the structure 
functions become scale invariant, they depend only on the dimensionless variable 
x B - 

The physical meaning of the structure functions is described by Feynman 
parton model which considers the nucleon seen from the scattering lepton as 
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Figure 1.5: Box diagram associated with the forward Scattering Compton amplitude 
in first order of a s . 



a gas of non-interacting, massless elementary particles, called partons. A field 
theoretic foundation of the parton model has been given by QCD thanks to 
its asymptotic freedom. At short distances and during short interaction times, 
QCD particles interact weakly and hence, if the nucleon is a bound state of rel- 
ativistic QCD particles, these appear in a DIS experiment as Feynman partons. 
Therefore the process should be described as a direct quark- virtual-photon in- 
teraction in the framework of QCD. Even if gluons should be considered partons 
together with quarks, they are not considered as partons in DIS since they do 
not directly interact with photon. However QCD Feynman model should present 
perturbative corrections, since quarks are not exactly free; as we mentioned in 
the previous section, perturbation theory induces logarithmic corrections to the 
pure scaling because of renormalization, but these corrections are computable 
in QCD. 

Neglecting QCD correction, one has the naive parton model in which: 

g a,P W a ,p{p,q) = - f <7 point (x B /y)qo(y) — , (1.20) 

where qo(y) is the parton distribution function 

<J P oint{xB/y) = e 2 S(x B /y - 1). (1.21) 
From this relation [27] one immediately has 

-W2(x B ) = =2F 1 (x B ) - e 2 q (x B ) 2 (1.22) 

x B 

The equality F 2 {xb) ~2xb F 1 (xb) 1 called Callan Gross equality, is experimen- 
tally verified. If one considers the radiative correction in QCD, and in particular 
diagrams with loops like the box (1.5), which has three fermionic lines and a 

2 e is the fractional electric charge of the considered quark 



1G 



gluonic one, one meets a divergence in the limit p 2 — > 0. This divergence arises 
in regions of momentum space, in which the gluon momentum is parallel to that 
of the massless quark. 

Considering all the first order in a s radiative corrections the cross section be- 
comes [27]: 

o point {z) = e 2 ■ [S(z - 1) + ^ (log(Q 2 /m 2 )P(z) + f(z))] , (1.23) 

where P(z) and f(z) are finite and log(Q 2 /m 2 ) contains a mass singularity in 
the limit m — > 0. 

The infrared problem is solved by the Altarelli-Parisi equations [6], which 
allow the inclusion of such a singularity in the effective density of the quark. 
Indeed the effective parton density depends by the scale //.: 

Qo(y) -> q{y, log(p 2 /m 2 )) = q (y) + Aq{y, log(p 2 /m 2 )), (1.24) 

where Aq(y, log(^ 2 /m 2 )) contains the mass singularity of the divergent diagram: 



...2, f , Mv) 

2tt 



>/- 

J x „ 



Aq(y, \og{^ 2 /m 2 )) = ^ log{n 2 /m 2 ) / dy^- ■ P(x B /y). (1.25) 



The parton density given in (1.24) satisfies the following evolution equation: 
d , a ol s {() f 1 q(y,t) 



dt 



q {xB ,t) = ^[ dy q -^-Pix B /y) + Oia s it) 2 ), (1.26) 
2tt J x - y 



whose solution depends on the initial condition at t = iqoix)). The trace of 
the hadronic tensor (1.20) in terms of the effective parton density becomes: 

W2(x B ,t) =-J^Jy^fl.e 2 [Six B /y-l) + ^itPix B /y) + .fix B /y))] = 
= -e 2 qix B ,t) +0(a s (t)) 

The function P(x) is known as the splitting function and it is evaluated sum- 
ming on the radiative corrections of the diagrams which contribute to the first 
order in a s of the cross section, isolating the collinear divergent part. We need 
a regularization in order to find the collinear divergent part of the diagrams. In 
literature the dimensional regularization is used and the collinear singularities 
appear as poles in e of mass-shell massless quark amplitudes in much the same 
way as UV-divergences do. In other words the dimensional regularization pa- 
rameter e has a universal meaning. In order not to mix divergences of different 
nature, the UV-divergences must be subtracted first. 

In our approach the quark is considered weakly bound inside the hadrons and 
hence off-shell. Let S 2 measure the off-shellncss of the parton initial state. This 
regularizes collinear divergences but not the infrared singularities associated 
with the single parton final states. This choice, which allows the absence of 
UV-divcrgence mixing, combined with the use of the Speer Sectors will allow 
us isolate the collinear divergences, which appear as divergent function of S 2 . 
The use of an off-shell amplitude together with the Speer sectors representation 
provides us an alternative method for the study of collinear singularities in QCD, 
which will be presented in this work. 

We consider this new approach worth a careful analysis, at least for physical 
and formal reasons, even if it is possibly more cumbersome than the standard 
general one. 
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Chapter 2 



Schwinger representation of 
Feynman amplitudes. 

In perturbative field theory the Feynman amplitudes are generally represented 
either in coordinate or in momentum space. 

There is another representation which is very important in practical cal- 
culation and will be systematically used in this work. This is the Schwinger 
parameter representation, that we shall introduce in this section. 

In the framework of this representation one can divide the parameter measure 
in Speer sectors and distinguish the singularities of the amplitudes more easily. 

2.1 Proper diagrams 

Scattering amplitudes and operator matrix elements within scattering amplitude 
states are related by L.S.Z. reduction formulae [28] to the Fourier transformed 
time-ordered Green functions, that is, to the Fourier transformed vacuum expec- 
tation values of time ordered products of fields and local composite operators, 
such as currents. 

In Feynman h~ expansion Green functions are related to series of con- 
nected Feynman diagrams. These correspond to sets of vertices connected 
by lines forming a connected polygon. 

Among the vertices one distinguishes those associated with the fields ap- 
pearing in the Green functions; in the diagrams they appear as end-points of 
single lines, the "legs" of the diagram. Computing physical quantities the lines 
connecting the field vertices to the rest of the diagram are omitted. The diagram 
is "amputated" of its legs. 

Considering a diagram contributing to a Fourier transformed amputated 
Green function, one has a product of terms some of which are propagators and 
some others appear as a momentum multi-integral of a product of vertex and line 
factors (propagators). Each of these terms corresponds to a subdiagram, which 
is one-particle irreducible (1-P.I.), also called proper, meaning that 
it is not only connected, but remains connected if one of its lines is omitted. 
The computation of the 1-P.I. diagrams is the technically difficult part of the 
calculations. Therefore we shall concentrate our study on 1-P.I. diagrams. 
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A 1-P.I. subdiagram has three kinds of vertices, all of them are end-points 
of at least two lines of the diagram. The vertices of the first kind correspond 
to local composite operators. Those of second kind correspond to Lagrangian 
interaction terms, which are connected through one line either to the other 
parts of the whole diagram or to a field vertex. The vertices of the third kind 
are only joined by lines of the same subdiagram. We call the vertices of the first 
two kinds external vertices since they carry momentum flowing into the 
diagram from outside. Internal vertices are those of the third kind, they 
do not carry any external momentum flow. 




Figure 2.1: Example of diagram with an internal vertex, two external vertices 
of the second kind, A and C, and an external of the first kind, B. 

If a connected diagram has / lines and V vertices, the number of loops L, 
i.e. meshes present in the diagram, is given by the following topological relation: 

L = I-V + 1. (2.1) 

A tree (sub)-diagram is a connected diagram without loop 

Defining P v as the total four-momentum flowing into the diagram through 
the vertex v, the four-momentum taken by internal lines will be determined 
identifying in the diagram L different loop subdiagrams, that is, closed not self 
intersecting polygons, and selecting a line in each loop (polygon). Notice that 
even if there arc lines which belong to more than one loop, one has to avoid 
selecting the same line for two or more loop subdiagrams or lines whose momenta 
differ for an external momentum. Asking for momentum conservation in all the 
vertices, the line momenta are determined in terms of the external momenta and 
of momenta carried by the loop lines, that we call loop momenta. A Feynman 
amplitude in d dimensional space-time is an integral of a product of vertex 
factors and propagators over the L loop independent momentum variables, hence 
it is an integral in L d variables. 

We start our analysis in a scalar theory, characterized by couplings with- 
out derivatives and scalar propagators, calculating the Schwinger representa- 
tion of 1-P.I. diagrams. Then we shall extend the results to the more general 
and physical framework of theories characterized by non-scalar propagators and 
couplings, hence to theories with spinor and vector fields and with vertices con- 
taining derivatives. 
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We limit ourselves to massless theories. 



2.2 Schwinger parameter representation of 
Feynman amplitudes in A0 4 theory. 

The scalar propagator corresponding to the z-th line, which carries a four- 
momentum ki in momentum representation is described by: 

(2.2) 



kf + ir] 



The scalar amplitude of a proper 1 diagram G with / lines, V vertices and 
L = I — V + 1 loops is the following: 



A G (P) = K (27r) d 5 d C£P)A G (P) 



= K{-i\) v ^ 



I r 



n 



(27r) d kf + it] 



n 

v=l 



(2.3) 

(2n) d S^(P v -J2^vik) , 



with: 



• d = 4 — 2e is the dimension of the space-time, with the e dimensional 
regularization complex parameter. 

• K is the symmetry factor of the diagram. 



lr The relation counts also for a connected diagram, but we choose to only deal with proper 
diagrams 
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• [i is the renormalization scale. 



• A is the coupling constant of the theory. 

• P = {Pi, P n } is the set of the external momenta which satisfy 



£)P„ = (2.4) 



1=1 

and P v is the total external momentum entering through the vertex v. 

• e v i is the V x I incidence matrix, where v refers to the vertex and I to the 
line: 

+1 v is the starting point of I 
e v i = { —1 v is the ending point of / (2.5) 
otherwise 

In order to pass from momentum to Schwinger parameter representation [5, 7, 
20] one has to make the following replacements: 



kf + irj 



dai e 



ia i (k l +27) ) 



(2Tr) d <5 (d) (p v - J2 tvlk^J = J d d y v exp -iy v (p v - ^ e vl h 



(2.6) 



Then for every loop one has a Gaussian integral in d-dimensions which is easily 
performed leading to an expression of the form: 



Aa(P)=^{-i\) vi — ZrW ! 

(4tt)— ,J; Uo 



X r r°° 1 exp[i(Q G (a,P)} 



don 



[P G (a)} d / 2 



(2.7) 



In this expression the Symanzik functions Pg{o) and Qg^tP) appear. 

Pg(ch) is a homogeneous polynomial of degree L in the a variables. It can 
be computed identifying the set Jt of all the tree subdiagrams of G which are 
maximal, that is, such that one cannot add a line to the subdiagram without 
closing a loop. This implies that all maximal tree subdiagrams are obtained 
deleting in G one line for every independent loop in all possible ways. Therefore 
the maximal tree subdiagrams of G, elements of Jt, have I — L lines. Notice 
that a maximal tree subdiagram contain all the external vertices of G. 

Once one has identified all the elements of Jt, one computes Pg(&) using 
the formula 

Tej T igT 

where the sum is over all maximal tree diagrams T contained in the set Jt- 
One also has: 

Q G {a,P) = -^ T -D G {a,P), (2.9) 

where Dc(a, P) is a quadratic form in the vertex momenta P v , whose coefficients 
are homogeneous polynomials of degree L + 1 in the a's. 

In order to select independent momentum variables, the set of the external 
vertices must exclude an arbitrarily chosen vertex, that we call vq. Indeed 
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P Va is not independent of the other P v 's since momentum conservation implies 

Therefore one has 

D G (a,P)= C G {a) vy/vo P v P v < (2.10) 

It remains to compute the coefficients Cc{a) v _ v i /v - 

For this one has to identify Jt 2 , the set of the 2-tree subdiagrams which 
can be obtained from the tree subdiagrams, deleting one of their lines. It follows 
that in general a 2-tree subdiagram is not connected and that, if it is connected, 
it does not contain all the vertices of G; in other words in general a 2-tree has 
two disconnected components and, more important, divides the set of vertices 
of G in two subsets of the vertices contained in each component. If the 2-tree 
has a single component, one of the sets contains a single element, the vertex not 
belonging to the 2-tree. 

Considering the set of the lines of G not belonging to a 2-tree, these can 
be identified with the cut lines when one cuts the diagram in two components 
breaking all its loops. We call this cut a complete cut. From this point of 
view the above mentioned two subsets of vertices are those lying on opposite 
sides of the complete cut. 2 




Figure 2.3: a) is an example of a 2-tree b) is an example of a tree c) is neither 
a tree nor a 2-tree 

Thus a 2-tree contains I — L — 1 lines and is the complement in G of a 
complete cut of L + 1 lines. Given two, possibly coinciding, external vertices w, 
v' and «o, one identifies a subset of 2-trees Jt 2 v,v' /v containing the 2-trees for 
which the pair v,v' and vq lie on opposite sides of the corresponding complete 
cut. 

One has: 

Ca(a)v,v/vo= E II a ' ( 2A1 ) 

Alternatively if we call C the set of complete cuts C, understood as a subset of 
lines of G, and C v y / VQ the set of complete cuts separating the pair v, v' from 
vq, one can write: 

C G (a) v,v' /vq 

cec vy/vo iec 

2 In the following with complete cut wc will refer to a set C of lines whose deletion breaks 
G into two tree subdiagrams each of which contains at least an external vertex. 
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The integral (2.7) might diverge when one or more parameters a vanish or 
tend to infinity. The first case corresponds to ultraviolet (UV) divergences; the 
second one to infrared (IR) divergences. 

In order to evaluate equation (2.7), one introduces the scale variable t ac- 
cording to: 

aj=*A- (2-13) 

It is apparent that the definition of the new parameters /? is not unique, 
since they are / coordinates of a J — 1-dimcnsional manifold which is identified 
with the quotient space of the positive sector of the /-dimensional Cartesian 
space and a positive semiline, whose points are identified by the scale factor t. 
This manifold is the positive sector of the I — T dimensional projective space. 

Therefore the integration measure becomes 

i 

Y[dai = dtt I ~ 1 diJ,(l3), (2.14) 
i=i 

where dfj,(p) depends on the choice of independent coordinates (3 of the pro- 
jective space and the polynomials Q , P e D depend on t in the following 
way: 



Q G (ai-t,P)=t-Q G ((3 l ,P), (2.15) 
P G {a l -t)=t L -Pg(Pi), 
D G (a l -t,P)=t L+1 -D G (p h P). 

In terms of the new parameters j3i, one has: 
AMP) U 2L H iX) v i—^ fdt M f e^tQ G {^P)] 

a G {p) - n ( i\) {47r)dL/2 j t t j nm [PGmd /2 

(2-d)L-i x v t T(I-dL/2) 



The argument of the Gamma-function is (-D/2), where D = dL — 21 is the 
diagram superficial divergence degree. In a four dimensional scalar theory D is 
an even integer. When D is a positive integer or null the function T(— D/2) 
has a singularity corresponding to the overall UV divergence of the diagram. 
It is important to notice that the amplitude can be divergent also when D is 
negative, since the diagram might contain divergent subdiagrams. 
As said above, after the integration over the scale parameter the Feynman inte- 
gral is transformed into an integral over a positive sector of (I — l)-dimensional 
projective space. One can choose among different possible parametrizations of 
the projective space. 

These are the most relevant for our purposes: 

1. the hypercubic parametrization, which satisfies the constraint J2i=i A = 1) 
whose integration measure is 

dm = (n^ ) 5 (!-I>); ( 2 - 17 ) 

\i=l J i=l 
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2. the parametrization based on Speer-Smirnov sectors, that is, on the divi- 
sion of the integration domain into subdomains. This will be described in 
the following. In each subdomain one of the (3 is assigned the value 1 and 
the others take positive values less than 1. 

The great advantage of Speer-Smirnov choice is that UV and the IR divergent 
parts of a diagram appear only in few sectors and hence if one is interested in 
singularities, one can limit the number of considered subsectors. 



2.3 Schwinger representation 
Feynman amplitudes. 



of tensorial 



If one considers a generic theory, like QCD, one deals with tensorial amplitudes 
in loop momenta because of the presence of fermion propagators, which differ 
from the scalar ones, since they have a line momentum dependent numerator fi, 
and derivative and hence momentum dependent interactions terms: the vertex 
with three gluons in QCD, shown in figure (1.1), introduces a function of the 
line momenta in the numerator of the amplitude. 

Therefore the generic amplitude in momentum representation, differently 
from the scalar amplitude, has as a numerator consisting in a sum of different 
rank tensors in loop momenta, which are contracted with tensors depending 
on covariant objects, like external momenta, gamma matrices and the metric 
tensor. Let's call the numerator N(P, k). 

The momentum representation amplitude for a proper diagram with / mass- 
less lines and V vertices is 3 : 



A G (P) = 



Wn[ 



(27r) d kf + it] 



N(P,k)]J 



v=l 



(2n) d S^(P v ~J2 e ^) 



i=l 



(2.18) 

In order to give a parametric representation of this integral, one has to proceed 
in the following way: 

• One chooses a line I for every loop of the diagram (the line which carries 



the independent loop momentum) and multiplies its propagator 



kf+iri 



by 



the factor e lkl ' Ul , where m is a vectorial parameter assigned to line I. 



One replaces N(P, k) with N(P, —i ■ d u ), which depends on the differential 
operators d u and acts on the scalar amplitude. 



The amplitude is [5]: 



^ L *N(p,-id u ) / J] 



1=1 

V 



iki-u t 



d d k. 



(27r) d kf + ir\ 



n 



(2TT) d S^(P v ~J2 e ^ 



\u=0 



(2.19) 



= ^(2ir) d 5 d (J2 P) N{-id u )A G (P, «)| u=0 



Cv is a factor according to Feynman rules 
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After performing the L gaussian integrals in d-dimensions, one obtains for 
Ag{P,u) the following relation: 

A G (P,u) = ^C v l ^l]J 
with 

Q G (a, P, u) = D G {a, P, u) = Q G (a h P„+V % v ^-)~ V ^. (2.21) 

P G (a) *Y 2ai ^i^i 

Equation (2.20), that one commonly finds in the literature, is not sufficient for 
our purposes. Indeed in order to profit completely of the Speer-Smirnov sector 
decomposition, it is convenient to interchange in equations (2.10) and (2.12), 
which defines the quadratic form D G , the sum over vertices with that over 
2-trees. 

Indeed, a given 2-tree T2, or better, the corresponding complete cut C, 
contributes in the vertex sum to the coefficients C G (a) v y / Vo , when v and v' lie 
on the opposite side to vq of the cut diagram; let us call Vt 2 , Vo the set of vertices 
opposite to vo- Then the contribution to D G from T 2 is equal to 

iec v,v'ev T2 ,v lec \ v 

It is therefore natural to consider the kinematic invariant variable (Y] v P v ) 2 = 
P^ 2 , where Pt 2 is the total momentum crossing the complete cut C. 
After this choice one has 

D G (a,P)= J2 P l ri a ' = E P cIl^' ( 2 - 22 ) 

T 2 EJt 2 l£T 2 cgc iec 

where Pt 2 = Pc is the total momentum crossing the complete cut C. 

Let us now consider D G (a, P, u), given in equation (2.21), and single out the 
contribution of a given 2-tree Tn and the corresponding cut C. This is given by: 

iec \ v v 1 ) 

In the sum over I', appearing in this equation, it is clear that the lines 
connecting pair of vertices on the same side of the cut C do not contribute, 
since the sign of t v y is opposite for the two vertices of the pair. Therefore the 
only lines contributing are the elements of the complete cut C whose end points 
belong to opposite connected components of G. These lines do not belong to T2 
and arc such that the union T2 U V is a tree subdiagram. Therefore the linear 
term in u in (2.21) is given by 

e e e ^ p ^ ■ ««' n a i ( 2 - 23 ) 

T 2 eJr 2 l'£T 2 1<£T 2 

T 2 ui'eJ T 



da. 



exp [iQ G (a,P,u )] , 
[Pg{u)Y^ 



|u=0> 



(2.20) 




2G 



ei' t T 2 is positive when u crosses C in the same direction of Pt 2 ; otherwise it has 
the opposite sign. 

This formula and the following ones can be written in terms of the complete 
cuts and the reduced cuts C, which are obtained omitting a line from 
complete cuts and are identified with the complements in G of the maximal 
trees. Therefore, together with the set of maximal trees Jt, we have the set of 
reduced cuts C, such that for any C E C one has aTgJr satisfying: 

C = G/T. 

Now the term in equation (2.23) is written: 

E ei',cPc-u v Y[ ai , (2.24) 

c/i'ec 

where we have identified ay ,t 2 with e;/ ; c, if C is the complement of T 2 . 

Now we come to the quadratic part in u of Da(a, P, u), for which the con- 
tribution of a given complete cut is 

Once again the lines contributing to this term are those belonging to C and 
connecting vertices lying in opposite sides of the cut diagram. Let us consider 
first of all the coefficients Uf -ui" for V ^ I". As above I' and I" are restricted by 
the same condition appearing in (2.24) and hence one has for the corresponding 
term in Do(a, P, u) in terms of the 2-trees and the complete cuts: 

\ XI € i',T 2 e;",T 2 ui, ■ ui„ J| ai = 

T 2 eJ T2 I'jti" i<tT 2 uiui' 

1'uT 2 EJt 
T 2 UI"£Jt 

^ E zv ,c zi" ,c up ■ up, Y[ a i- ( 2 - 25 ) 

cec i'jti" iec/(iui') 
c/i'ec 
c/i"ec 

Considering instead the terms in uf, one has two contributions of opposite signs 
in (2.21), the quadratic part in Qq and the last term. The first contribution is 
positive and, following the analysis of the other terms, turn out to be equal to 

i e e n ( 2 - 2g ) 

T 2 EJt 2 l'£T 2 i£T 2 UZ' 

i'ur 2 eJ T 

The second contribution is 

-jftwEj-iEn-Ej- 

i' Tej T leT v 
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Combining these terms we find 

-iE n 



Oil. 



(2.28) 



TeJ T l'<£T 1<£TUI' 



Indeed the first term given in (2.26) selects the maximal trees which coincide 
with T 2 U I' for some T 2 . These are all the maximal trees containing I'. Since 
the expression given in (2.27) refers to all the elements of Jr, the difference 
corresponds to maximal trees which don't contain 

The formula in terms of the reduced and complete cuts becomes 

-i ee-' n ««■ ( 2 - 2 9) 

ceci'ec i^i'Aec 

In conclusion we have this expression in terms of maximal trees and 2-trees: 

D G (a, p, u )= p T 2 n ai + e e ^ pt * ■ ^ n ai + ( 2 - 3 °) 

T 2 eJr 2 1<£T 2 T 2 eJr 2 ,i, 1<£T 2 

(T 2 uheJ T ) 

( \ 



e u * n ai ~ e Ui ■ u j e ^ n ai 

i<£T 2 uhuij 



T£j T ,i l<£TUh T 2 Gj T2 ,i^j 
(h<£T) (l t UT 2 eJ T ) 



) 



While in terms of the complete and the reduced cuts one has the alternative 
formulation: 

D G {a,P,u) = J2 P cI[ a i+ E ^■,cPc-u l \{a l + (2.31) 

cec lec cec,i, lec 
(hec) 
(c/hec) 

( \ 



E u i 1 1 ai ~ E Ui ■ u i ei > ce i< c H ai 

cec,i iec/h cec, iec/(hui 3 ) 



(hec) 



cec,i^j 
{c/hec) 
(c/i 3 ec) 



I 



We recall the equivalent form for (2.22), 



cec lec 



the other Symanzik function being: 



p G(a) = EIl a <- 

dec led 



(2.32) 



(2.33) 



Going on we shall use this alternative form. 
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Considering equation (2.18), one can write a generic amplitude as a lin- 
ear combination of homogeneous polynomials of the components of the exter- 
nal momenta whose coefficient are given by suitable integrals in the Schwinger 
parametric space. Indeed let us decompose N(P, k) into a sum of products of 
homogeneous terms in the external momenta P and internal loop momenta k: 

N(P,k) = J2K(P)-K(k). (2.34) 

n 

Let djyi and be the mass dimensions of the factors. 

Correspondingly in equation (2.20) one has N(P,—id u ) = X] n ^n(^) ' 
N^(—id u ). The action of N^(—id u ) on cxp i Dc p ^'^ \ u=0 generates a sum 

of contributions corresponding to the possible actions of the u-derivative either 
on the linear or on the quadratic part in u of Dc(a, P, u). The number of times 
the w-derivatives in iV* act on the u-quadratic part of Dc(a, P, u) must be an 
even number, 2a, since the result is computed at u = 0. Therefore acts on 
dfji — 2o times on the n-linear part of Dc(a, P, u). 

One has correspondingly a factor p ^|"| a with R{a) homogeneous of degree 
(L — I) a from the quadratic part of Dc(a, P, u) and a factor — s ^ °d P ^--i a with 

Pa(a) N " 

S(a, P) homogeneous of degree d N ^ — 2a in the momentum components and of 
{d N ^ — 2a)L in a from the linear part. Therefore one has a global factor of the 
type: 

M 

where we denote by [X] the integer part of the positive number X, and 
P^ a (a,P) is a homogeneous polynomial of degree d^i — 2a in the external 
momentum components and of degree 

(d N , n - 2a)L + (L — l)a = L d N * n — (L + l)a (2.36) 

in a. 

The a-degree is apparently non-negative. Therefore P^ a {a, P) is a linear 
combination of a-monomials with momentum-dependent coefficient: 

PU^ P ) = J2 C Ur(P)tl aXl ' T - ( 2 - 37 ) 

T 1 = 1 

Consequently in the parametric integral the numerator N, whose dimension dpj 
is directly related to the global degree of the amplitude G, 

d G = d N + dL-2I, (2.38) 

corresponds to the factor : 

""^ nf-, (dual'-) 

fc=0 a=0 t r G\^) 
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where we have denoted by b the degree in k of the homogeneous polynomi- 
als N^(k) in the decomposition (2.34). Furthermore Ol 6 jT is a homogeneous 
polynomial in the momentum components of degree 



de = djv — 2a 



and one has, as shown above, 



^ Xi tT = L(b — a) — a. 



(2.40) 



In a scalar theory, (b = 0, a = 0), this quantity is null as expected. 
At this point in much the same way as in scalar theory, one introduces the 
scale variable t according to (2.13) and integrates on it. 

One obtains the analogous generic relation to the scalar case, (2.16): 



■I+L-Ld 



^f- d N -2a 



where 



^ > t a=0 6=0 



I G ,(a,b,T)(P) = » ^(^a) 



and 



Ld 

tJ-r,a = I ~ -y - a - 



(2.41) 



(2.42) 



(2.43) 



In chapter 4 we shall study the ultraviolet and infrared divergences related 
to lG,(a.b,r){P) decomposing the integral into the Speer sectors. 

2.3.1 An example: the vertex correction 




k in 



One computes the off-shell, dimensionally regularized amplitude of the diagram 
in figure using the expressions (2.20) and (2.31). The amplitude in momentum 
representation is 



J (27T 



d d k J p @+fay"fa P 



fa" 4+ fa" % 
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(2.44) 
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where the numerator depending on k is due to the two fcrmionic internal lines, 
and one has Cy = —ie 2 s ecp • , since there are two QCD vertices gluon-fermion- 
fcrmion and a QED one. 

One assigns the fourvector u to line 1 and from (2.8) and (2.31) one has: 



Pg{o) = Qi + a 2 + a 3 ; 
r> < d ^ D G (a,P,u) 



(2.45) 



Pa{a) 



1 / u 2 

n -, s \D G (a, P) + a 3 p ■ u - ai q ■ u - 

Pg(oc) V 4 



with D(j(a, P) — p 2 aia 3 + q 2 a x a2 + So^o^ and S = (q + p) 2 . 

Passing from momentum to parameter representation, the amplitude be- 
comes: 



AJf = -/i 2e CV 



2i 



( 47r )d/2 



(P G (a)) d / 2 



( a 3 ^ - a 2 ^ ) 7 P 



exp [iQ G (a, P, uj\ 

[P G (a)] d/2 = ° ~ 
(qi + a 3 ) + i> a 3 ) ry'' 



P G (a) 2 



Pc(a) 



At this point one scales the parameters with respect to the variable t, like in 
(2.13), and, after i-integration, one obtains 



2Cv r 



(Pg(PW 

P 2 PiPz+q 2 Pifr + SfcPz 



(ft J*-/3a ^)7 M W (ft + ft) + * ft) 



(p 2 ftft +q 2 ftft + S /3 2 ft) 



(2.46) 



Here the second term is apparently UV-divergent since it diverges in the limit 
e^O. 

The MS renormalized amplitude is given by 



7 M C v 
e (47r) ; 



(2.47) 



where we have subtracted the UV-divergence due to the pole of T(e) — \ + C- 
The result is: 



" ~~ 73^ Jo 
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(2.48) 



- 7 M log 



(p 2 ftft +g 2 ftft+Sftft) 
P G (ftV 



P 2 ftft + g 2 ftft + Sftft 



We will study the infrared divergences of this amplitude once introduced the 
Speer sector parametrization. 
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Chapter 3 



Singularity families and 
Speer sectors 



In my degree thesis [1 7] I have studied the scalar Feynman amplitudes in the space 
of Schwinger parameters, divided into Speer sectors, showing that this division 
separates the amplitudes in the sum of terms with different singular behaviour. I 
have identified and classified the infrared divergences for scalar amplitudes. The 
construction of the Speer sectors was proposed by Speer in the article [8] in 1975. 
In this chapter we present the construction of the Speer Sectors due to Smirnov 
in [10]. The sectors one obtains are slightly different from the those identified 
by Speer, but related to them. They are more numerous than the Speer ones and 
the range values of the parameters are limited all between and 1, instead of 
taking different range values. 

The advantage of Smirnov 's choice is the availability of power counting for- 
mulae 

3.1 Preliminary definitions 

• An irreducible diagram is a diagram which either is 1-P.L, but it is 
not the union of two 1-P.I. diagrams with a single vertex in common, or 
it consists of a single line. 



Figure 3.1: a) is irreducible, b) is 1-P.L, but not irreducible. 

• The parts of a diagram are the maximal irreducible subgraphs of the 
original diagram. 




a) 



b) 
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• A link [8] in a diagram in a massless theory is a connected diagram 
which contains all vertices of the diagram and is minimal: any subgraph 
obtained after the deletion of any part of the diagram is not anymore a 
link. 




Figure 3.2: a) is a link b) is a link c) is not a link since the subdiagram is not 
minimal 



3.2 Construction of a singularity-family 

Let G be an irreducible Feynman diagram with / internal lines and L loops. A 
(singularity)-family T in G is a set of subgraphs of G, which are either 
links or irreducible subgraph of G, and which do not overlap, that is, either 
they have no lines in common or are contained into one another. This set of 
subgraphs is also called a forest. 

Since the diagram G is a link, we include it as the first element of the family. 
Then we delete any line l\ in G. Let us denote by G/h the subdiagram of G, 
obtained deleting the line l\ from G. Gjl\ contains all the external vertices of 
G, however in general it is not a link since it is not necessarily minimal; G/h is 
the union of parts, some of which form a link. This link, that we call the second 
link of the forest and we denote by (L ,r coinciding with G) and the other 
parts are new elements of the singularity-family. 

Then we delete a further line li in L\^j=. If some parts of Li^/1% form a link, 
L>2,f, we add L,2,r and the other irreducible parts of L\^jli to T and continue 
deleting a further line, I3, in L^.?. We repeat the same procedure until the last 
link, Lk,r is broken by a further deletion. Let us notice that each of the first 
k deletions has opened a loop in the corresponding link whereas the k + 1-th 
deletion has broken into two disconnected subdiagrams, each containing 
external vertices. At this stage the family is composed of k + 1 links contained 
into one another and a set of trivial or non-trivial irreducible subdiagrams which 
do not intersect, that is, which have no lines in common. A trivial subgraph 
contains a single line. We choose an arbitrary non-trivial irreducible element of 
J 7 and delete one of its lines, Z 7 , adding the parts of j/lj to T . This procedure 
is repeated considering at each step the set of non-trivial irreducible elements of 

which are minimal, i.e., which do not contain other elements of J 7 , choosing 
an arbitrary element of this set and deleting a line in it. 

The procedure stops when one is left with only trivial elements. Now the 
construction of the family is complete. Since the deletions of the first set have 
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opened k loops contained in the links L,^, whereas each deletion of the second 
set has opened a loop in a minimal non-trivial irreducible element of T and the 
sequence of deletions stops when there are no unbroken loops left, the number 
of deletions is equal to L + 1. 

Furthermore in every non-trivial element of T we have selected a line. The 
set of selected lines corresponds to a complete cut. It is apparent from the 
construction procedure that the unselected lines belong to the trivial elements 
of the family and are identified with them. Therefore the total number of trivial 
and non-trivial elements of T is /, the number of lines in G. There is a one-to- 
one correspondence between the elements of JF and the lines in G, l-y,F- 

What is worth noticing is the fact that the sequence of deleted lines identifies 
a partially ordered complete cut in the diagram and hence the set of undeleted 
lines identifies a 2-tree. Indeed the first set of deletions including the one 
which breaks the last link, the k + 1-th deletion, forms a cut which separates 
the set of external vertices in two subsets. Let's call it the breaking cut 
Cjr, corresponding to the family J- . This is completely ordered. The number 
of lines contained in the breaking cut equals the number of the links in T . 
However the cut is not complete; indeed non-trivial irreducible subdiagrams 
remain after the considered deletions. 

We show an example of construction of a family for the five loop diagram of the 
figure: 




If one deletes in order the lines 1, 2 and 3 in T, the first element, 
three links are formed: L(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), L(3, 4, 5, 6, 7, 8, 9, 10, 11, 12), 
and 1/(4,5,6,7,8,9,10,11,12). Deleting the line 4, there are not links 
anymore; so one includes in the family the parts of the subdiagram: 
#(7), #(10), if (5, 11,12), and #(6,8,9). Deleting the lines 5 and 6 in the not 
trivial parts, one finally has #(8), H(9), #(11), #(12) and the procedure stops with a 
family of 12 elements: 

r e L(r/{i» e L(r/{i, 2}) e L(r/{i, 2, 3}) e #(7) e #(10) © (3.1) 

ffi#(5, 11, 12) ® #(6, 8, 9) #(8) © #(9) © #(11) © #(12) 

The 2-tree associated with this family is formed by the lines of the trivial parts of 
the families, that is, {7, 8, 9, 10, 11, 12} and the breaking cut is {1, 2, 3, 4}. 

In summary we have met three types of cuts: 

• The complete cut, that is the set of L + 1 not ordered deleted lines, is 
the complement in a diagram G of a 2-tree; each of two tree subdiagrams 
contains at least an external vertex. 
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• The reduced cut is the complement of a maximal tree in a diagram G and 
contains L not ordered lines. 

• The breaking cut, corresponding to a family J 7 , is an ordered set of lines, 
whose number depends on the singularity family. 

Given a family T, generated as described above, one has a corresponding 
complete cut Cjf, which is identified with the set of deleted lines. 

We want to show that, given any family T and any complete cut C, C 
intersects all the non-trivial elements of J 7 ; more precisely, for every 7, non- 
trivial element of T, there is a distinct line in C, which is contained in 7, 
meaning that one can set a one-to-one correspondence between the L + 1 lines 
in C and the non-trivial elements of J 7 . 

In order to prove this, let us notice that, given the set of L + 1 non-trivial 
elements 7 of J 7 , it is possible to identify, in general not uniquely, L independent 
loops in G, asking that every loop is contained in a 7. This is possible since the 
elements 7 do not overlap. The complete cut C contains at least a line in every 
loop, otherwise it should not be a complete cut; however there is a one more line 
in C and one more non-trivial element in J 7 than there are loops in G. One has 
to select a line lk in C which is contained in the smaller link Lk and such that 
C /Ik is a reduced cut. This is always possible at least in one way. Indeed, first, 
the cut C intersects Lk since it separates the external vertices, while Lk joins 
them. Second, at least for one line lk in the intersection of C with Lk, C/lk is a 
reduced cut. Indeed T2 i c ; the complement of C , separates the external vertices 
while Lk joins them: this means that there is at least one line in C fl Lk, which 
added to T 2 gives a T e Jt, since this line joins two disjoint subgraphs of T 2 . 

Once the reduced cut C/lk is extracted from C and lk is associated with Lk, 
which contains it, C/lk cuts all the loops and hence it contains one distinct line 
for every loop. Since there is a correspondence between loops and non-trivial 
elements of J 7 containing them, we have identified for every non-trivial element 
of J 7 a distinct line of C contained in it. 
An example can make this argument clearer. 

Let's consider the diagram G of figure (2.1) and two different sets of non-trivial 
elements contained in two singularity-families, T a and Tb, related to the diagram G. 

Deleting the lines {4, 6, 2, 3} respectively one obtains the four links shown in figure 
(3.3). 




Figure 3.3: These are, together with the diagram itself, the non-trivial elements of a family J-" a , 
produced by the first L — 4 choices. In this case they are all links. The remaining two elements of 
the family, which are trivial, are obtained by deletion of any line of the last link. 
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We see that each of the above deletions produces a link. So one has as many links 
as there are loops, if we don't consider G. In order to complete the family one has 
to choose a line between the ones of the last link L 4 , for example line 5, so that the 
diagram breaks into two trivial parts. In this example the breaking cut coincides with 
the complete cut C? a : {4, 6, 2, 3, 5} 

Alternatively one can delete from the diagram G the lines {4,5,3}. One has the 
non-trivial elements in figure (3.4), showing links together with irreducible diagrams. 




Figure 3.4: These are the non-trivial elements of the family J-fj. The deletion of the lines {4, 5} 
produces two links. From the second deletion one has also the part H(6, 7). From the choice of the 
line 3 one has another irreducible clement, H(l. 2). 

The breaking cut contains the lines {4, 5, 3}, since the choice of line 3 breaks the 
last link. In order to complete the family one has to delete one of the lines {1,2} of 
H(l,2) and one of the lines in H(6,7). 

Now we verify that it is possible to create a correspondence (not unique) between L 
loops and L non-trivial elements of a family, by selecting in every non-trivial element 
an independent loop. Indeed for the family T a one can establish for example this 
correspondence 

G o {2, 4, 5}; LI O {6, 7}; L a 2 o {1, 2}; L% O {3, 5, 7} (3.2) 

and for the family Th'- 

G -H- {1,4,5}; L,\ o {3, 5, 7}; H(l, 2) o {1, 2}; H (6, 7) o {6, 7}. (3.3) 

In this way we have created a correspondence between 4 loops and 4 non-trivial ele- 
ments of the two families. 

Now we consider a complete cut C = {1,3,4,5,7}. It contains 5 > 4 lines. As 
explained above, we first construct from C two reduced cuts related to the two families. 
One has to take the last link of the family and consider the lines I, which are contained 
in the intersection of the last link and the cut. For the two families one has 



L% nC= {1,5,7} 
L|nC = {l,3}. 



o ^ _ n (3-4) 



In case a one has to exclude both lines 1 and 7 since neither {3, 4, 5, 7} nor {1, 3, 4, 5} 
are reduced cuts; their complements contain loops. Hence one must choose l a = 5. In 
case b one must choose l a = 3. 

Therefore, given C, in this example one has only one reduced cut in both families: 



C a = cjja = {1,3,4, 7} 
C 6 = C/? = {1,4,5,7}. 



6 _ nnb m , = » ( 3 - 5 ) 



Now, taking into account of (3.2) and (3.3), we construct one of the possible corre- 
spondences between the lines of the complete cut and the non-trivial elements of the 
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family: 

G o 4; I?«T; L a 2 o 1; L§ ^> 3; L;«5 

G <-> 4; L^o 5; L§ -H- 3; #(1, 2) 1; H(6, 7) o- 7, 1 ' 

where we have associated i| = 5 with the last link of the family T a and l 2 = 3 with 
the last link of the family Tb- This shows how we create a correspondence between 
distinct lines of a complete cut and non-trivial elements of a family. 



3.3 Speer-Smirnov sectors 

The crucial step in our program is the identification of the decomposition of 
the positive sector of the projective parameter space seen in the previous chap- 
ter into subsectors, each associated with a particular parametrization of the 
projective space. 

This is accomplished referring to the singularity families T . We have seen 
that a singularity family is a partially ordered set of non-overlapping subdia- 
grams of G including G itself and that there is a corresponding complete cut 
Cjr which singles out a line in each not trivial element of the family. 

The partial order of the elements of J- is determined by their inclusion prop- 
erties: every subdiagram precedes in the order the elements of T contained in 
it. With each element 7 of J 7 we have associated a unique line l 7 ,jr which is 
contained in 7 and in no element of T contained in 7. If 7 is non-trivial this 
is a line of the complete cut associated with the family, Cjf. The set of the 
lines associated with the links, which are elements of the breaking cut, have a 
natural complete order since the links are contained into one another. The lines 
associated with other irreducible elements are partially ordered, as already said, 
and they are ordered after those of the link. 

Therefore the partial order in T corresponds to a partial order in the com- 
plete cut, which consists in the ordered lines of the links and in the partial 
ordered lines of the non-trivial irreducible elements. Since every line of G is 
associated with a particular clement of J 7 , the family contains as many elements 
as there are lines in G, as we have said. 

Notice that the lines of the 2-tree, which is the complement in G of Cjf, are 
elements of the trivial subdiagrams in T . Having recalled these points, we come 
to the parametrization associated with a family T. One associates a parameter 
tj, which takes values between and 1, with every element 7 di J 7 , while one 
sets to equal to 1. 

The parameter /?/, appearing in the scalar amplitude (2.16) and in the generic 
one (2.41), is set equal to the product of the t 1 corresponding to the elements 
of T containing the line I. That is one sets: 

a= n *t ( 3 - ? ) 

The range of values taken by the /3's after this prescription defines a Speer 
sector. It is apparent that, due to the range < t 7 < 1, the partial ordering 
of the elements of F corresponds to a partial ordering of the /3's. In a given 
non-trivial element of the family T the line with the largest /3i parameter is the 
line Z 7i jf. 
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Coming back to the first example, we see that the Speer sector corresponding to 
the family (3.2) is: 

(e*5 > an, ai2 
CK6 > C*8, ag (3.8) 
ai7, am 

The biggest parameter is assigned to line 1: ai = tr = 1. Calling i; the variable 
corresponding to the i-th element of the family, starting from L(F/{1}), the parameters 
have the following explicit form: 

ai = tr = 1, OL2 — ti, as = titz, e*4 = £1*2^3, «5 = t\titzt§, e*6 = t\titzt-r, 

07 = tit2t'it4, aio — tit2tgtr > , an = i^isWio, <*12 = tit2t3tetii! (3-9) 

C*8 = il^2^3^7^8, = tltzt^tv- 

It is not difficult to verify that different sectors do not intersect, since dif- 
ferent families correspond to different sectors. Indeed considering the family 
construction procedure, there is necessarily an irreducible diagram, e.g. G it- 
self, with which different families associate different lines. This corresponds to 
a different ordering of the lines of the irreducible diagram. 

Furthermore one verifies that the union of all the sectors covers the whole 
positive sector of the projective space. 

In order to show this second point, let us consider any ordered choice of the 
Schwinger parameters > on 2 > ... > aj 7 > 0. Let us identify with the 
scale parameter t, and hence choose: 

ft, =l>ft 2 >... >&,. (3.10) 

We have to show that this corresponds to a Speer-Smirnov sector. 

From the construction procedure it is clear that k ± is Iq ■ The situation with 
li 2 is less simple; indeed U 2 could either belong to a link contained in G/l^ or 
to a non-trivial irreducible part of G/l± or else to a trivial irreducible part. In 
any case l i2 is identified with the line of either the link or the irreducible part. 

The above chosen order is compatible with this identification. Let us now 
consider ?j 3 ; once again it must be contained in one of the minimal elements 
introduced into the family after deletion of lines and l i2 , let k 3 be the line 
associated with this element; let us furthermore reduce what remains of this 
element after deletion of k 3 in its irreducible parts, we can now consider Z,- 4 
repeating the same procedure. Once reached kj, we have identified a family 
■^ii,i2,..,i/j an d hence a Speer-Smirnov sector, compatible with the considered 
order of parameters. This proves that the sector decomposition exhausts the 
whole parameter space. 

In order to make the construction clear we discuss a particular example. 
One can then verify our claims concerning the disjointness and completeness of 
sector decomposition. 

Example of construction of the complete set of singularity- 
families for a two loop diagram 

Let's consider the two loop diagram of the figure with three external vertices 
and an internal one and let's build the entire set of singularity families following 
the procedure above. 
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Figure 3.5: Scalar two loop diagram with three external vertices and an internal 
one 

From the initial diagram every choice of the first line produces a link after 
its deletion. So the links which are formed are so many as the the lines of the 
diagram, which are five: 

L(l,m,i,j), L(m,n,i,j), L(l,m,n,j), L(l,m,n,i) and L(l,n, 

Let's start considering the first link L(l,m,i,j). In order to obtain another 
link, one has to delete anyone of the lines except for j; in this way the links 
L(m,i,j), L(l,m,j) or L(l,i,j) are formed. The construction stops choosing 
anyone of the lines of the link and one obtains different families with the same 
structure which we call of type "A". If otherwise one chooses line j in the 
initial link what remains is a non-trivial irreducible subgraph, H(l,m,i), which 
has to be included into the family. The procedure stops choosing one of the 
lines contained in H(l,m,i). These families have a different structure, since 
they contain a not banal irreducible element. We call them of type "B". 

Thanks to the symmetry of the diagram the families which contain the sec- 
ond link are obtained from the above described ones exchanging the lines i < — > j 
and I < — > n . 

From the third link L(l, m, n,j) only the deletion of line n produces families 
of type A, since the link L(l, m, j) is formed. Otherwise from the deletion of line 
I, which produces H(m,n,j), one has three families of type B. Finally deleting 
a further line, either j or m, one adds to the family the new link L(l, n) and a 
further trivial element, either m or j. The procedure stops choosing a line in 
L(l,n). We call this set of families type "C". We distinguish this set from set 
A since its elements contain a two line link. 

The families which contain the fourth link are obtained for symmetry from 
the ones which contain the third one. 

Finally from the fifth link one can alternatively delete line l/n or line i/j. In 
the first case the link L(n,i, j) /L(l,i, j) is formed and so one has families of the 
type A; in the second we have families of type C, since the deletion produces 
L(l,n)®j/i. 

In the end the number of the families is 54: 30 of type A, 12 of type B and 
12 of type C. 
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They are shown in the table (3.1). 



Table 3.1: Singularity families 
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Table 3.2: All the singularity families related to the diagram in figure (3.5) are 
shown. Every family contains also the diagram itself, which in this table is 
omitted. 
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Chapter 4 



The Symanzik functions in 
the Speer-Smirnov sectors 

In this chapter we will show that the Speer-Smirnov sector division plays an 
important role in the study of the asymptotic behaviour of the amplitudes. In 
every sector a power counting is possible for the study of the infrared divergences. 
This is shown for a generic theory. 

In order to evaluate exactly the divergent part of the amplitude, we will use 
the Mellin-Barnes transform, which will be described in the following chapter. 



4.1 Symanzik functions in a Speer-Smirnov sec- 
tor 

We have seen in chaper 2 the expression (2.41) of the amplitude Aq{P) in 
Schwinger parameter representation, given by a sum over the indices a, b, t. 
We want to focus on a single term of the sum, in particular we recall 



and express it in the Speer-Smirnov parameter representation. 

In this representation Ia.(a,b.T){P) decomposes into the sum of different con- 
tributions Iq ( a b T )(-P), each of which is evaluated in a Speer-Smirnov sector: 

lG,(aAT)(P) = Y, I G,(»AT)( P ) ( 4 - 2 ) 

s 

with 

rS (p\-u 2L 'Y(J — a) f dfXS nf=iA A< ' r u *\ 



,(a,6,r) \ r ) ~ t 1 L \ 1 2 > J Z^W-I+b n I- 

In the last expression the measure d/is and the Symanzik functions Ps and D$ 
are labelled by S, since they arc different in different sectors. 
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From (3.7), specifying the sector parametrization, we can deduce the corre- 
sponding integration measure of the sector amplitude, which is 

i 

= n = n ^ (7)_i dt ^ ( 4 - 4 ) 

1=2 7 G^ 

where the product runs on all the elements of the family, except for the diagram 
G itself, and 1(7) is the number of lines contained in the element 7 of the family. 

From the relation which defines the Symanzik functions (2.33) and (2.32), 
changing the variables from (3 to t 7 according to the definition of the sector, one 
has the general expressions for the Symanzik functions in a sector: 

p s - EIK (7n(5) ; 

°s = E P cIK (7nC) ' ( 4 - 5 ) 

where we recall that C is a complete cut, C is the reduced cut, which contains 
L lines (one line less than the complete cut), I(j n C) is the number of the 
lines in the intersection between an element of the family 7 and a cut C and 
7(7 n C) is the number of the lines in the intersection between an element of 
the family 7 and a complete cut C. 



In the previous chapter we have seen that a singularity family J 7 , and there- 
fore a Speer-Smirnov sector, identifies a particular complete cut Cjr. It is associ- 
ated with the kinematic invariant Pq f , which is the square momentum crossing 
the cut. 

Therefore it follows that 



D S = II 4 (7n ^ } [ p c, + E P c II *r 




where 



M7)=-f(7nC)-/( 7 nC f ). (4.7) 



It is important to notice that 60,^(7) has integer non- negative values and the 
condition 5c,r(j) = for every 7 6 T identifies C with Cjr. 

In order to obtain a similar relation for Pg , we have to introduce the ordered 
reduced cut Cjr. Cjr is obtained omitting from Cjr the line Ijr associated with 
the smallest link in T . 

The expression for Ps is: 

p s =u t^ ncv) ( 1 + e n 4 e ^ (7) ) = n < (7nc " ) n *(*7) (4.8) 

where 

j e/ ( 7 )=^nC)-/( 7 nC r ), (4.9) 
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which has non-negative integer values and it is identically null in F if and only 

if C = Cjr. 

At this point let's come back to the relation (4.3) and express it in the new 
variables t 1 . Combining the above equations together we have: 

ILJP) = ^ m - f - «) / II d h ih_ 7 J+6 (4-10) 

with 

f (7)-,6,r = Hi) - Hi n Cr)(b -a+l) + I{l* k n 7 )(i| - I + a) + £ i£7 A,. T . 

In the last expression I(ljr Cij) — I(-f n Cjr) — I{"f n Cjr) can either be 1, if Ijr 
is in 7 or null if 7 does not contain Ijr. In other words it depends on whether 7 
is a link or not. 

It is important to notice at this point that our formula (4.10) holds true 
for dimensionally regularized unsubtracted amplitudes. Concerning the UV- 
subtracted amplitudes, if the UV- divergence is primitive, there appear loga- 
rithms together with rational integrands. The case of amplitudes with UV- 
divergent subdiagrams in the sectors where these divergences appear is discussed 
by Smirnov in [10]. 

We shall not consider anymore UV-divergent amplitudes since we shall dis- 
cuss in this thesis a one loop example in which the !7V-divergences are primitive 
whenever they appear and hence do not give any contribution to the collinear 
divergent amplitudes. 

Therefore forgetting UV-divergences we notice that, if one considers the 
IR-safe theories mentioned in the introduction, that is, theories in which the 
interaction corresponds to operators with mass dimension strictly equal to 4, 
such as QCD, computing the amplitudes with non-exceptional Euclidean exter- 
nal momenta, one has absolutely convergent integrals which are analytic in the 
kinematic invariants. Indeed this follows from the known IR-power counting 
theorems [8, 13, 29]. 

This result is very important for us since it guarantees that the exponents 
£(7)0, b,r in equation (4.10) are positive and hence the measure 

TT £ £ (7)a,f,,T-l 

df a , b . T = LI dU f 7 J +1 ) 777 ( 4 - 12 ) 

is positive and integrable in the I — 1 dimensional hypercube < i 7 < 1. 
Therefore we have for the sector amplitude the expression: 

I S a,,r(P)= [dfil b / {I ~t~ a) - (4-13) 

J A 1 ^-* 

From equations (4.13) and (4.6) we see that T (P) is regular and analytic in 
the Euclidean region if Pq < 0, and hence A5 < 0. A5 vanishes in boundary 
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points of the / — 1-dimensional hypercube together with Pq ■ These are points 
where some t 7 vanish. This is the starting point of our singularity study. 

However, to complete this study, we have to change our point of view. In- 
stead of looking at a given sector- (family), we must consider one, or more than 
one, kinematic invariants whose vanishing together characterizes the limit under 
study. We have seen that once the sector- (family) is given, this identifies a com- 
plete cut and hence a kinematic invariant coinciding with the square momentum 
crossing the cut. 

If we start from a set of kinematic invariants with the same (or possibly 
proportional) negative value — £ 2 , we have to identify the set of families J^, 
the corresponding sectors and complete cuts Cjm , such that the kinematic 
invariants P^ ? belong to the above set and hence are equal (proportional) to 

— £ 2 . We call Q the set of the complete cuts associated with kinematic invariants 
vanishing together with — £ 2 . 

With these definitions the singular part of the amplitude appears in the sum 
over the sectors of T (P) : 

iUt(p) = e / d A,r < » N tr,°L-^- a ' (4 - w) 

where 

N si {t,) = l+ II 4 C '" 5(7) (4-15) 
cec 6 

is a positive polynomial and 

M s( (t,)=J2 P c II 4°" f(7) ( 4 - 16 ) 
c<£c ( 7 e^« 

is a negative polynomial which in general vanishes in boundary points of the 
hypercube. 

Our purpose is to analyze the behaviour of j| b T when ^ 2 — > and this can 
be done employing the Mellin-Barnes transform. 

However one can study the limit of the amplitude in every sector before using 
the Mellin Barnes to verify the presence of the divergence sector by sector. On 
this purpose we study the example of the two loop diagram, whose singularity 
families have been derived in the previous chapter. 

4.1.1 An example 

Suppose we are interested in the limit k 2 — > of the amplitude corresponding 
to the diagram (3.5). We consider the other kinematic invariants of the same 
order of magnitude: q 2 = p 2 = u 2 . The four complete cuts, through which k 
flows, are: 

{n,j,i}, {n,j,m}, {n,j,l}, {n,m,i} 

and are drawn in figure (4.1.1). 

We show that the amplitudes associated with the corresponding families- 
sectors might be singular in the k 2 — > limit. 
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Figure 4.1: Complete cuts through which k flows 

Notice that each cut corresponds to many families. We select few examples. 
Let's start considering the type A family: 

Ta,a = G ® L(l,m,i,j) © L(l,i,j) ®l®i, (4.17) 

shown in table (3.1) whose corresponding complete cut Cjr Al is {m,n,j}. We 
assign tc = 1 to the element G and the variables t*, with i = 1,2,3,4 to the 
elements of the family in the order they are written. 

Since we are only interested in the possibly singular behaviour of the am- 
plitude in the k 2 — s> limit we replace II5 defined in equation (4.12) and N$ 
defined in equation (4.15) by 1. 

Concerning MgA,i(tj), defined in equation (4.16), we know that its dominant 
terms are related to the complete cuts C, for which 5c = X) 7 ^ M = 
1(7 n C) — 1(7 n Cjr A 4 ) is minimal, for example the cut {to, n, /} into which p 
flows for which Sq = 1. 

We stop looking for other complete cuts, since for all of them 5c > 1 and 
hence they correspond to negligible corrections. Thus the contribution of this 
sector to the amplitude is given by: 



/ 



t 2 dt 2 dt< 



(4.18) 



k 2 + t 3 cu 2 ' 

Apparently it diverges in the limit k 2 — > 0. 

Another family of type A, which is worth analyzing , is 

Ta,\ = G ® L(l,m,i,j) ® L(m,i,j) ® m ® i. (4.19) 

corresponding to the complete cut Cj? A 1 — {n, I, j}. 

The dominant terms in M S i,i(t 7 ) correspond to the contributions from the 
cut {n, i, j}, through which p flows, and the cut {n, m, j}, through which q flows. 

Therefore the amplitude is approximated by: 



/ 



t 2 dt 2 dt 3 dU ^ 2Q > ) 



k 2 + {t 3 + U)uj 2 ' 
which doesn't diverge. 

As a last example let's consider a family of type B: 

F B ,8 = G® L(l, to, n, i) © H(l, m,i)®l®i (4.21) 
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The complete cut associated with this family is C? B 8 = {j, n, m} and the 
dominant terms of M S B,s(t 7 ) correspond to the complete cut {j, l,m}, through 
which p flows, and to {j, I, i}, through which q flows. 
The scalar amplitude is approximated by: 



and in t± diverge. 

At the end of the analysis one finds that the scalar amplitude diverges in 
the limit k 2 — > in 20 sectors. 

An analogous conclusion can be reached considering the limit q 2 — > 0. 

On the contrary in the limit p 2 — > the only families which potentially give 
divergent contributions are of type A, and are associated with the complete 
cut {I, m, n}. One can verify that there is no collinear divergence in the scalar 
amplitude in p 2 — > 0. 

At this point we need a mathematical tool which makes us compute exactly 
the divergent part of the amplitudes, without approximations. 




(4.22) 



and has a stronger divergence in the limit k 2 



— > 0, since both the integral in ti 
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Chapter 5 



Mellin-Barnes 
representation 

The Mellin-Barnes representation [15] allows the analytical evaluation of the 
divergent infrared part of the amplitude, since it separates the singularity in the 
kinematic invariant from the rest of the amplitude. In this chapter we consider 
this important mathematical tool and in the next chapter we will apply it to a 
physical problem. The method we want to show is based on Speer-Smirnov sector 
division combined with the Mellin-Barnes transform. It is a successful method 
in the analysis of collinear divergent amplitude. 

5.1 Definition 



The Mellin-Barnes transform is a useful mathematical tool based on the follow- 
ing formula: 
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where a and b are real non- vanishing numbers with the same sign, a; is a non- 
vanishing complex number and the path V [—zoo, +ioo] is a continuous line 
going from Zi = —iT to Zf = iT in the T — > oo limit, leaving the poles of r(cr), 
z n = —n with n G (0, oo), on its left-hand side and the poles of T(x — a), 
z m = x + m with me (0, .., oo), on its right-hand side. 

This formula is easily proved if b/a ^ 1 since the integration path can be 
closed, without further contributions, enclosing all the poles of T(cr) if b/a > 1 
and all the poles of T(x — a) if b/a < 1. 

The formula can also be generalized to the case of complex non- vanishing a 
and b, with a + b / 0, turning the integration path to V [—e 1 ^ oo, +e 4 ^ oo] , with 
< 4> < arg (^) ± | < n and the same conditions concerning the poles of the 
integrand. 

5.2 Employment of the MB transform in a Speer 
Sector 

We shall exploit (5.1) in order to evaluate the behaviour of an Euclidean 
Feynman amplitude given in equation (4.14) when £ 2 — > 0. 

Formally, choosing a = ^ 2 A^ S {(i 7 ) and b = — M s «(t 7 ), we have: 

= A £ /<U (-A/ S <(i,)r»-' (5-2) 

JV[-ioo,+ioo] 1 \ ? ly SZ\h)J 

We have seen how the integration path should be chosen depending on the value 
of the ratio (~ g2jv^(t\ ) • However, in order to profit of the Mellin-Barncs, we 
have to change the order of the integrals considering 

lUr - I ^-(e)-°r(a)T(I-^-a-a) (5.3) 

JV[-ioo,+ioo] Z7rt L 

Taking into account that N S ((t 7 ) is positive and bounded in the hypercube 
(0 < t 7 < 1), if — A'/ S {(t 7 ) were bounded from below by a positive number, 
for small enough £ 2 the two formulae in equation (5.3) and in equation (5.4) 
would coincide. This is however not the case, since M s s(t 7 ) might vanish in the 
boundary. 

A consequence of this vanishing is the fact that the integral 

m - J m$,,r (-Ms^)r^- 1 (-|^f£f . (5.4) 

instead of being an analytic function of cr, has a finite number of poles. 

It follows that, while in the analytic case the path in equation (5.4) would 
be closed around the poles of r(cr), in the second case the path must leave not 
only the poles oiT(a), but also those of /(cr), on its left-hand side and those of 
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r(J — — a — a) on its right-hand side and the amplitude must be computed 
closing the path around the poles of T(cr) and those of /(cr). 
Let us make this more clear presenting a trivial example: 



The pole from the rr-intcgral transforms that of T{a) in a — from first to 
second order and hence produces the term — log(£ 2 ) which gives the singular 
part of the original integral. The remaining power series gives log(l + £ 2 ) which 
is the analytic part of the integral. 

We shall give a less trivial example of the use of (5.1) applying it to a DIS 
process at one loop. 




(5.5) 
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Chapter 6 



Collinear singularities in 
DIS 



In the first chapter we discussed the important physical process of DIS. Here we 
present an off-shell strategy for evaluating analytically the collinear IR-divergent 
part of the structure functions of DIS to first order in the strong constant as, us- 
ing the parametric representation of the Speer sectors combined with the Mellin- 
Barnes transform. 

6.1 The computing strategy 

In the framework of the naive parton model the hadronic tensor is computed 
in terms of the parton (quark) distribution function qo(y) and of the partonic 
tensor W^ u (yp,q) using the formula: 

W^(p,q)= [ dyq {y)W%{yp,q), (6.1) 
Jo 

where y is the fraction of proton momentum taken by the quark. 

In QCD the naive partonic tensor is computed from the tree approximation 
forward quark-virtual-photon Compton scattering by means of: 

Tr^^+<j)j u } \ 
(p + q)*+iV ) [ ' 

Thus the partonic tensor is related to the absorptive part of the tree approx- 
imation forward Compton amplitude in figure (6.1), according to the optical 
theorem. 

Taking into account the radiative corrections, one has to consider the forward 
virtual photon-quark scattering amplitude corresponding to the diagrams in 
figure (6.2). 

In all these diagrams the initial state and the final state coincide and consist 
of a quark with momentum p and of a virtual photon with momentum q. 

A scattering amplitude depends on six independent kinematic invariants. 
These are the square momenta (masses) of the four external legs, that is, of 
initial and final, possibly virtual, particles, the square center of mass energy, the 



W^(p,q) = —Im 



1 

4tt' 
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Figure 6.1: The forward tree- approximation Compton amplitude 




Figure 6.2: The first order in as forward amplitudes in DIS. 

Mandelstam variable S, and the square momentum transfer, the Mandelstam 
variable T. 

We denote the independent kinematic invariants of these amplitudes by p 2 , 
the square momentum of the quark, q 2 , the square momentum of the virtual 
photon, S = (q+ p) 2 and of course T = 0. 

Considering the off-shell quark-virtual photon forward (Compton) ampli- 
tude, we have negative virtual photon square momentum (q 2 ), positive center 
of mass energy S and null T, while we choose negative virtual quark square 
momentum (p 2 ). 

In principle one should foresee a singularity in the amplitude at T = 0. How- 
ever this is not the case for helicity /angular momentum reasons. A vanishing 
angular momentum in the crossed (T) channel q + q — >• 7* + 7* is excluded by 
helicity conservation and consequently the amplitude is regular at T = 0. It is 
not analytic for S > 0, as expected, since we are studying the absorptive part 
of the amplitude. 

However if we complete our off-shell choice considering an unphysical neg- 
ative S together with the negative square momenta of the virtual scattering 
particles, we expect analyticity for the general reason mentioned in chapter 4, 
that is, infrared power counting. Therefore the general idea is the following. 

• We start computing the amplitudes with negative p 2 , q 2 and S, where the 
amplitudes are analytic. 

• We separate the singular part when p 2 — > 0~. This is expected to diverge 
as log(— p 2 ) and in fact it does. We compute the coefficient of log(— p 2 ) 
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in the Mellin-Barncs expansion of the amplitude. This gives the complete 
p 2 -expansion of the amplitude, however we disregard the regular contri- 
butions in the present analysis. 

• This coefficient of log(— p 2 ) is expected to correspond to an analytic func- 
tion of S, in the region of negative S and q 2 . By the Mellin-Barnes formula 
and sector decomposition we compute explicitly the analytic coefficient of 
log(— p 2 ) (as a matter of fact it is the sum of analytic contributions) and 
we analytically continue it in the Re(S) > region, where we find the ex- 
pected branch-cut. We compute the discontinuity, which is directly related 
to the singular part of the trace of the partonic tensor under study. 

We repeat this analysis for each graph, noticing however that, on account of 
L.S.Z. formula [28], diagrams d and e have to be divided by 2. 

For each spinorial amplitude we select the mass-shell divergent term con- 
tributing to the trace of the hadronic tensor and we compute the trace of the 
spinorial matrix multiplied by^, since we want to sum over the helicities in the 
(logarithmically singular) mass-shell limit. 



The box 




k 1 n 



n3 k 




k-p 

a 



In momentum representation and in the Feynman gauge the amplitude is 



e 2 e 2 s cf 
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r d 4 i 
J W 



d 4 k TrVYM+jthM 



Y ([k 2 ) 2 ( q + k) 2 {k-p) 2 y 



(6.3) 



and we have 



W^\p,q)=Im{A^{p,q)). 



Let's start contracting the amplitude with the metric tensor in order to single 
out the collinear divergent box contribution to the trace (in Lorentz indices n 
and v) of the partonic tensor. This gives us the radiative corrections to the 
linear combination of structure functions in the first line of equation (1.18). 
The radiative corrections to the second line are given contracting the expression 
in equation (6.3) with p^p v . However we shall see that this does not contain 
singular contributions in p 2 — \ 0~. 
Computing the trace we get: 



TrtyYftYU+^lvftlp] = 16 (2k ■ pk ■ q + k 2 (k ■ p - p ■ q)) . 



(6.5) 
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At this point we pass to the Schwinger parametric form following the proce- 
dure shown in chapter 2 and disregarding the terms proportional to p 2 in the 
numerator, which give vanishing contribution in the p 2 — > limit. 
We get 

with K = as e l CF . 

Then, integrating over the scale factor t and hence passing from the a to the 
ftparameters, we get: 

,(«) v f MP) (, 2ft \ 1 

with J = , fi ^ 

1 -D a (AP) = P 2 ft(ft +ft) + g 2 ft(ft +ft) + ftftft + ir?. { ' 

We have shown explicitly the infinitesimal imaginary term ir], which accounts 
for the time-ordering in the Feynman amplitudes. The amplitude is the sum of 
different terms corresponding to different Speer sectors. We want to single out 
the terms which diverge in the limit p 2 — > 0. To this purpose we have already 
said that we can focus on a subset of sectors in which the scalar amplitude of 
the box is singular. This amplitude is proportional to 



A (a) 



DC 



DaW,P) 2 

and decomposes into 12 Speer sectors corresponding to the parametrizations: 

a) ft = sf> sft ft = s a, ft = sf3 <H- s, ft = 1; 

b) ft = 1 sft fe = sa,ft = s^ 1, ft = s; 

c) ft = 1 s, ft = s a, ft = s <-> 1, ft = sft 

d) ft = 1 sft /? 2 = s, ft = s,3 1, ft = sa; (6.9) 

e) ft = s o sft ft = 1, ft = s/3 <-> s, ft = sa; 
/) ft = sft ft = 1, ft = «*, ft = s; 

g) ft = sft ft = s, ft = sa, ft = 1, 

where each of the first five lines corresponds to two sectors, giving the same 
contribution, because of the symmetry for exchange of the parameters ft and 
ft • 

It is easy to verify that the box scalar amplitude is divergent in the limit 
p 2 — >• only in the six sectors of a) b) and c): 



n\ f 1 da ds ( f 1 da.\ f 1 _ 

a > JO (p 2 +a(q 2 s+S)) 2 ~^ \h a 2 J JO (q 

h \ r 1 da ds p 2 ->0 / r-l da \ rl ds 

U ) JO (p 2 +a(q 2 + S s)) 2 ~^ \J0 TP J Jo (q 2 + Ss) 

„\ fl da ds d£ / r-\ fa\ r-l 

C > Jo ( P 2 f}+a(q 2 +S sf3)) 2 ~^ \J0 a 2 J JO 



S) 2 ' 



( 9 2 +Ss/3) 2 • 

Notice that, following the method discussed at the end of section 4.1 in order 
to identify the families and hence the sectors accounting for the singularities in 
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the limit of vanishing p 2 , that is, selecting the families for which the momentum 
crossing the complete cut Cjr coincides with p, we should have selected sectors 
a) and b) forgetting sectors c). 

However we are not facing a counter-example to the above method; indeed 
we have to consider, together with the families a) and b) associated with p 2 , 
the families c), which correspond to the kinematic invariant T, since T is null: 
= T < \p 2 \. 

Now, once selected the interesting sectors, we consider again the expression 
for the amplitude given in (6.7), which is the sum of two terms. The second 
term, which is proportional to /?2, does not contribute to the collinear divergence, 
since in all the considered sectors (3 2 — sa vanishes with D a ((3, P) for p 2 = 0. 

Thus the singularity is only due to the first term in the integrand of (6.7) 
and hence we must evaluate the sum of the contributions to 

_ P q J Pa{P) 2 D a {^P) 

from the six sectors a), b) and c), that is 

sdsdadB ( 1 



A {a ^-2Kva( sdsdadp f 

~ P "Jo (l + s(l + a + P)) 2 { 



(1 + /3)(p 2 +q 2 sa) + Sa + irj 



1 1 

+ 



(1 + j3 s)(p 2 + q 2 a) + S sa + it] (1 + s)(p 2 (3 + q 2 a) + S s a /3 + if] 

In order to single out the collinear divergence we apply the Mcllin-Barncs for- 
mula, introduced in chapter 5: 

A^ = -2Kp-q f £(V)-n^r(l-a) (6 . 10 ) 

JV[-ioo,+ioo] Zm a 

sdsdpSa ( (1 + 13)-* 



o (1 + s(l + a + f3)) 2 \[-a(q 2 (l + l3) S + S)-i V ] 1 - a 
{1 + 13 3)-° (l + s)-' 7 



[-a(q 2 (l +j3s) + Ss)- if]} 1 '' 7 [-a(q 2 (l + s) + S s p) - irj} 1 ^ J ' 
Let us perform first the a integral; this has the form: 

/ —a a A(a), (6.11) 
Jo « 

where A(a) = ^ +s ^^ +a ^i is analytic and non-vanishing in the integration 

domain. It follows that the above integral is +TZ(a), where TZ(a) is analytic 
in the whole complex cr-plane. 

This pole in a — is the singularity considered in chapter 5 which must 
be enclosed in the path encircling the T(a) poles. Thus the cr-integrand has a 
double pole in a = 0, while this pole would have been simple in the absence of 
collinear singularities. 

Notice that if we had used dimensional regularization as IR-regularization, 
the new pole would have been in a = e. In this case, performing the cr-integral 
one would obtain a term — - from the T-pole in a = and ^ — from the new 
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pole. Then one could "go to the limit" p 2 — > 0. Choosing, arbitrarily, e real and 
negative, the second term would vanish in the limit, leaving a \ singularity, as 
a memory of the collinear one. 

We believe that our procedure is mathematically more sound. 

Thus, applying the above considerations to equation (6.11), computing the 
cr-integral and selecting the contribution of the double pole, the rest giving 
analytical contribution in p 2 w 0, we get: 



4? } 



2p-q 2 f 1 sdsdp 

K - log(-p z ) 



div -q 2 ° ' Jo a + s(l+f3)) 2 \(l+(3)s + ^-iV 



+ ^ + 5 . (6.12) 

l + Ps+fs-irj l + s + ps(3-iri) 

This is an analytic function of the variable with a branch-cut on the negative 
real axis. We compute the discontinuity and we set 

4 = 1--, (6-13) 

which now is a negative quantity, since S is positive 1 , and we also set z = 
l + s(l + <* + /?). 

The discontinuity is: 



W% [a) 



r 1 r 1+2s dz 

= ttK \og(-p 2 ) ds — (6.14) 

div J J 1+s Z 

[8(xz- 1) + 5(xz- s) + 5(-z(l -x) + l + s)]. 



Notice that S = q 2 + 2p ■ q = q 2 - 3- and that = i. Hence ^xi = l. 
After the integration we have 

f «) {^)e(\-x)6{x-\) + {^)e(x-\)e{i-x) 

K \og{-p 2 )ir- I b) x6{x)6(\-x) +(l-2x)6(+\-x)8(x-\) 
( c) (l^)9(x)9(+l-x), 

(6.15) 

where we see that different sectors correspond to different physical regions. 
Finally the sum is 



W P,{a) 



2,(1-*) 



= KTrlog(-p 2 )^— -L (6.16) 

div Z 



which is the box contribution to the collinear divergent imaginary part of the 
amplitude, we were looking for. 



Before studying the other diagrams, we still have to contract the box ampli- 
tude (6.3) with p»p v . 

The numerator in momentum representation is proportional to: 

Tr^Y¥M+¥Wl P ] = -32(k-p) 2 (k-p + p-q). (6.17) 

1- We set x = xb, Bjorken variable. 
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One has 

. -2K iP . 9 fj MI3) m±^±M (6 . 18) 

It is very easy to verify that the amplitude is not divergent in the limit p 2 — > 0. 
Indeed, in analogy with the above case we have a /3|-proportional term. /3f 
vanishes with D a (j3,P) 2 in the (p 2 — > 0)-limit. 

The vertex correction 




Let's recall equation (2.49), in which we evaluated the first order in a s amplitude 
of the vertex correction: 



ecp roc 

2vr Jo 



(Pom 3 



7 M log 



(p 2 PiPz + q 2 /3 1 p 2 +S + irf) 



+ 



y p 2 Pife+q 2 Pifc + Sfcfo, 
We can write the amplitude of the diagram b in terms of A^: 

e Tr$Y¥+4)^) 



Afj v {p,q) = -i 



4tt 



S + iri 



(6.19) 



In much the same way as for the box amplitude, we first contract with the 
metric tensor and then with p^p v . 

Contracting equation (6.19) with the metric tensor, multiplying by two, since 
the diagrams in the figure give the same contribution, and finally disregarding 
the terms which are proportional to p 2 , we have: 



2A 



(*>) 



K- 



2p-q J <W) 



S+ir, J P b (f3P 



D b {P,P) 



log 



D b (J3,P) 



with 



D b (fi,P) = P 2 l3 1 l3 3 + q 2 f3 1 f3 2 + Sf3 2 f3 3 + i V . 



(6.20) 
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There are 6 Speer sectors corresponding to the following parametrizations: 



a) 


Pi 


= s,(3 2 = 


sa, p 3 


= i; 


b) 


Pi 


= 1,P2 = 


sa, Pz 


= s; 


c) 


Pi 


= sa, p 2 


= S,P 3 


= i; 


d) 


Pi 


= sa, Pi 


= 1,03 


= s; 


e) 


Pi 


= S,P 2 = 


1,& = 


s a; 


f) 


Pi 


= 1, P2 = 


S, P 3 = 


s a. 



(6.21) 



It is easy to verify that the scalar amplitude, corresponding to diagrams b and 

is UV-fmite. It has IR-divergences in the limit p 2 — > corresponding to the con- 
tributions from the sectors a) and b), since in this limit Db((3, P) is proportional 
to a. Therefore we have to compute (6.20) in the sectors a) and b). 

In this equation we have put into evidence two terms. The second one 
comes from the subtraction of the UV-divergence. It is collinear finite, since it 
is a parametric integral of the logarithm of a function which vanishes on the 
boundary of the integration domain in the p 2 — > 0-limit. 

In the first term we pick out the only contribution which is collinear divergent 
in the sectors a) and b), which is: 

2At = K —Jim^ 1+M ^MP7n- (6 ' 22) 

The rest is finite since it is proportional to p2, which vanishes together with 
Db(P, P) in much the same way as for the box. Summing the contributions 
from the two sectors, we get: 



2A^ = -K^(S- q 2 )!im^W 



-p 2 + a(—q 2 s — S) — ir\ 
—p z + a(-q z — b s) — IT] 



We apply the Mellin-Barnes transform and we consider the collinear divergent 

part. 

This is: 



(*>) 



div 

( \ 



if] 1 + Jr s - if] f 



(6.24) 



Again this is an analytic function of the variable with a branch-cut in the 
negative real axis and a pole in — irj. 

The novelty of the vertex correction contribution lies in the presence of the 
pole superimposed on the branch cut. Thus computing the absorptive part 
we should take into account the branch-cut discontinuity and the contribution 
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proportional to —iirS(S) coming from the pole. This corresponds to a single 
scattered parton final state which is also present in the collinear-divergent con- 
tributions from diagrams d and e. 

However the term iw5(S) multiplying the vertex correction in diagrams b 
and c is ill-defined, since we see from equation (6.24) that the coefficient of the 
Dirac delta contains the ill-defined integral L ds/[(l + s) 2 (s — i-q)] . 

A general remark is here in order. Even if the Fourier transformed Feynman 
amplitude corresponding to a given diagram is the product of the contributions 
from the parts of the diagram, this product structure is not suitable for spectral 
analyses. Indeed for this purpose one has to introduce a spectral representation 
for the whole amplitude and compute, e.g. its imaginary part, on the basis of 
this representation. 

In the light of this comment it is easy to verify that the two above mentioned 
inconsistencies compensate each other. Indeed let us look more carefully at 
equation (6.24). It can be written as follows: 



2AP 



= K 



log(~p 2 ) 



div 



o s (1 + sy \ s + 4 -it] jk-iv, 



+ 



S/q 



i r 1 sds i r 1 

'--ivJo (1 + sf jis-irj Jo 



ds 



(6.25) 



Now, considering the imaginary part and setting 

^ = 1-1 

q z x 

we have: 



2W; 



P,(b) 



div 



= Ktt log(— p 2 ) x 



+ 



f dt (r~ 

J 1/2 \l-X 



1/2 



dt 
l-t 



t(5(x-t)-6(x-l)) 



-S(x - (1 - t)) + S(x -t) + 6{x - (1 - t)) 



-S(x - 1) log(2) 



where we have changed the integration variable t = 

The question is if this is a distribution in < x < 1. To check this point, 
we multiply 2W^'^ \ div by a C^ 1 j function cj){x) and integrate over x. 
We get: 



Jo 



dx 2W; 



P;(b) 



div 



Ktt log(— p 2 ) 

.1/2 



dx ■ 



1/2 



l-x 



Mx) - 0(1)) 



da; _£_0( X )_ (105(2-1)^(1) 



j\x X -^ 



+ 0(1). 



(6.26) 
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Defining as in literature the distribution jj^j^ as 



Jo (1-*)+ Jo 1- 



x 



we identify equation (6.26) with: 



Kir log(-p 2 ) 

diu \ (1 — X) + 



+ S(l-x) 



We sum this result with (6.16) getting 



(W£ (a) +2W£ (b) ) 



From (1.20) and (1-23) one computes the splitting function P(x) 



P(x) = c F [ -1+JL- + 25(1 -x) + C 6(1 -x] 



(1-x). 



(6.27) 



(6.28) 



(6.29) 



(6.30) 



where C is the coefficient of term from the diagrams d and e. Notice that 
diagram / does not contribute to the collinear-divergent part of g a, ^W a ^ since 
it is proportional to the tree diagram term multiplied by log(5). 




q.u 



Considering that, on account of L.S.Z formula, diagrams d and e have to be 
divided by 2, one has 



1 as e 2 CF 



log(— p 2 ) xS(l — x). 



div 4 7T 

Hence C = — ^ and we get the known value of the splitting function: 
P (a0 = CF (^ + §*(!-*)). 



(6.31) 



(6.32) 



We have still to contract 2A$ V with p^p v . One easily sees that the result 
vanishes in the collinear limit. 
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Conclusions 



We have seen in the last chapter how one can study collinear divergences in DIS 
cross sections by considering the quark mass-shell limit of Feynman amplitudes 
and exploiting the simplifying power of the Speer-Smirnov sector decomposition 
of the parametric integral expressions for the amplitudes together with that of 
the Mcllin-Barnes technique. 

In particular, we have shown that the off-shell regularization is a natural 
tool for the analysis of collinear divergences in inclusive cross sections, and that 
the use of this regularization is strongly simplified by the Speer-Smirnov sector 
decomposition. The recourse to Mellin-Barncs formula is the obvious final tool 
for the computation of singular parts of the amplitudes. 

The suggestion of initial state off-shell regularization of inclusive cross sec- 
tions as a reasonable "physical option", as an alternative to the widely adopted 
dimensional scheme, is one of the results of the present thesis. However the 
principal goal of this thesis was to identify a general method of analysis of mass 
singularities in the Feynman amplitudes of massless field theories, which could 
be automatically applied to multi-loop and multi-external-vertex diagrams. For 
this reason we have described with great care the role of different kinds of cuts 
in the construction of the parametric form of a Feynman integral. These cuts 
are sets of lines in Feynman diagrams whose deletion either breaks the orig- 
inal diagram in two parts, or transforms it into a tree diagram, or else does 
both things together (chapter 2). We have also analyzed the role of the cuts 
in the identification of the Speer-Smirnov singularity families (chapter 3) and 
sectors (chapter 4). The idea was that the identification of sectors contributing 
to the IR-singular parts associated with the vanishing of kinematic invariants 
is a central step in the study of these singularities, and that the method could 
be extended to multi-loop and multi-external- vertex diagrams developing suit- 
able software tools. As a matter of fact, in the last year an example of such 
tools (Fiesta [30]) was developed by Smirnov and Tentyukov. However, to our 
knowledge, the method has not yet been used in the explicit analysis of IR 
divergences. 

A systematic use of the techniques presented in this thesis may in principle 
prove useful in different contexts of phenomcnological interest. High-order per- 
turbative calculations in QCD, for example, may take advantage of a systematic 
separation of singular contributions. From a different point of view, an efficient 
way of identifying IR divergent terms to all orders in perturbation theory can be 
useful when one is faced with the problem of resumming the whole perturbative 
series in special kinematic regimes, where powers of large logarithms of ratios 
of invariants spoil the reliability of fixed-order calculations. 

There are, of course, many points that need a further and deeper discussion. 
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For example the gauge independence of the singular parts of the forward am- 
plitudes [14], has been mentioned in short occasional remarks. This could be 
better justified extending the basic argument proving gauge independence on 
the charged particle mass shell of massive QED amplitudes [31]. Indeed in any 
gauge theory there is a general connection between the gauge independence of 
the "physical" amplitudes and Ward-Slavnov- Taylor identities [32]. In the case 
of amplitudes involving only local physical (BRS invariant) operators one finds 
directly that their partial derivative with respect to any gauge fixing parameter 
vanishes. The case of amplitudes involving charged fields, which are not BRS 
invariant, is more difficult, indeed one finds that the same derivatives are less 
singular on the mass-shell than the original amplitudes; poles are replaced by 
branch cuts. The idea which must be verified in a general situation is that in 
mass-less theories the same thing happens since any gauge parameter derivative 
of a collinear divergent amplitude is less singular, or even null in the mass-shell 
limit. This is fairly clear in our DIS example, indeed, using Slavnov- Taylor 
identities, one shows that the Feynman-gauge parameter derivative of the non- 
amputated forward scattering amplitude is proportional to the box diagram 
amplitude in which the gluon propagator has been replaced by a scalar propa- 
gator. The result is apparently proportional to the box amplitude (there is a —2 
factor), therefore it is proportional to log(— p 2 ). However to get the derivative 
of the forward scattering amplitude one has to amputate the amplitude, which 
means multiply it by p 2 , and hence the collinear divergence disappears and the 
gauge derivative vanishes in the mass-shell limit. 

A second point that we have not discussed with sufficient completeness is 
the effect of UV-divcrgent parts on the IR analysis. This is essentially the 
consequence of the choice of the Smirnov version of Speer's construction, which, 
due to the crucial role given to the links, particular subdiagrams containing all 
the external vertices, is more suited for the study of IR divergences. 

In the original Speer construction IR and UV divergent subdiagrams partic- 
ipated in the singularity families on the same footing. A further study on this 
point should recover this advantage of Speer's construction which is crucial if one 
tries to work with high-order diagrams. As a matter of fact one should compare 
Speer-Smirnov singularity families with Hepp-Zimmermann forests and show 
that the UV divergent parts correspond to singularity families and hence sec- 
tors containing one-particle-irreducible UV divergent components consistently 
with UV power counting. Therefore, being clear that Speer's sector decom- 
position is useful in the singularity analysis, while it does not seem particu- 
larly suited for regular diagram calculations for which other sector decompo- 
sition have been developed together with a suitable software, we believe that 
our work has been convincing enough concerning the advantages of the Speer- 
Smirnov sector decomposition combined with the Mellin-Barnes formula in the 
IR-singularity analysis, also in view of the study of multi-loop multi-leg am- 
plitude. This point is also stressed by A.V. Smirnov, V.A. Smirnov and M. 
Tentyukov in few recent works [33]. 
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